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I.  INTRODUCTION 


Applied  Research  Laboratories,  The  University  of  Texas  at  Austin 
(ARL.UT),  has  been  conducting  research  under  Contracts  N0001 4-82-K-0425 
and  N00014-86-K-0176  to  investigate  both  theoreticaliy  and  experimentally  the 
basic  properties  of  the  optoacoustic  sound  source.  This  report  pertains  to  the 
latter  contract. 

Several  mechanisms  of  sound  generation  may  be  identified  under  the 
general  title  of  optoacoustic  sound  sources.  The  simplest  one  is  the 
thermoacoustic  process  which  employs  direct  heating  of  the  acoustic  medium  to 
produce  a  controlled  local  thermal  expansion  which  in  turn  generates  sound 
waves.  1  The  thermal  energy  is  delivered  by  a  laser  beam  to  the  water  surface 
without  mechanical  contact.  A  moving  thermoacoustic  source  may  be  obtained 
by  deflecting  the  laser  beam  with  a  rotating  mirror,  causing  the  position  of  the 
thermal  sound  source  to  change  as  a  function  of  time.  This  allows  the  signal  to 
be  Doppler  shifted  much  more  than  it  can  be  with  conventional  sound 
projectors,  which  have  to  be  in  physical  contact  with  the  acoustic  medium.  The 
study  was  pursued  from  two  directions,  the  fundamental  physical  processes  and 
the  feasibility  for  sonar  applications.  The  results  of  the  investigations  have  been 
report  ed.2.3.4 


Although  it  had  been  speculated  that  the  moving  thermoacoustic  source 
may  be  more  efficient  than  the  stationary  source,  the  feasibility  study  found  the 
efficiency  upperbounds  to  be  identical  for  both  cases.  The  upperbounds  were 
found  to  be  constrained  by  the  physical  properties  of  the  medium  and  the 
energy  output  of  the  thermal  source.  Furthermore,  the  conversion  efficiency 
upperbounds  were  found  to  be  too  low  to  be  useful  at  sonar  frequencies  at  the 
levels  of  laser  energy  that  are  considered  to  be  practical. 3  To  obtain  higher 
efficiencies  without  requiring  impracticaliy  high  levels  of  laser  energy,  the  focus 
of  the  study  has  shifted  under  Contract  N00014-86-K-0176  to  nonlinear 
optoacoustic  processes.  Nonlinear  processes  have  been  found  experimentally 
to  be  more  efficient, ^  but  they  are  not  well  understood. 


1 


There  have  been  a  number  of  publications  on  this  topic.  A  large  proportion 
of  the  reported  works  have  been  for  cases  where  most  of  the  energy  is  spread 
over  several  megahertz.®-^  The  results  reported  by  Hickman®  were  in  the  range 
of  sonar  frequencies  but  they  were  very  sparse,  and  the  experimental 
conditions  appeared  to  be  imprecise.  The  most  extensive  set  of  results  were 
reported  by  Maccabee  and  Bell.®'^®-^^'^^  jhejr  data  covered  a  wide  range  of 
operating  conditions.  Their  results  showed  significant  finite  amplitude  distortion 
and  extra  attenuation  of  the  acoustic  pressure.  They  used  a  "rocket"  model  to 
estimate  the  reactive  pressure  on  the  water  surface,  and  linear  superposition  to 
construct  the  acoustic  pressure  signal;  they  have  found  that  the  model 
predictions  are  not  consistent  with  the  experimental  results. More  suitable 
blast  models  are  those  of  A.  N.  Pirri^^-^^  who  successfully  used  them  to  model 
the  momentum  transfer  between  a  high  power  laser  blast  on  the  surface  of  a 
solid.  Unlike  some  models  that  have  been  reported  that  are  specific  to  solid 
surfaces,^®  the  Pirri  models  are  general  enough  to  be  adapted  for  laser  induced 
blasts  over  water. 

The  theory  of  the  nonlinear  optoacoustic  generation  process  has  only 
been  partially  analyzed  in  existing  publications.  Wu^^  made  a  theoretical 
analysis  using  linear  superposition  and  a  model  similar  to  Pirri's  model  but 
without  any  experimental  verification.  Lyamshev  and  Naugol'nykh^^  produced 
a  similar,  but  one-dimensional,  model.  Feiock  and  Goodwin^®  also  produced  a 
one-dimensional  model  which  included  a  numerical  solution  of  the  pressure 
wave.  In  all  cases,  they  ignored  backscatter  and  absorption  of  the  laser  energy 
by  the  water  vapor,  and  the  effects  of  the  deformation  of  the  water  surface. 

Sigrist  and  KneubuhP®  made  excellent  experimental  measurements  of  the 
axial  acoustic  pressure  in  the  nearfield  but  did  not  proceed  beyond  the  linear 
thermoacoustic  model  of  Hu.20  Their  equipment  did  not  allow  the  measurement 
of  beam  patterns,  and  therefore  it  was  not  possible  to  analyze  the  spatial 
character  of  the  sound  source.  Maccabee  and  Bell®'^®  have  made  extensive 
experimental  studies  but  they  have  not  gone  beyond  the  "rocket"  model  and 
linear  superposition  in  their  analysis.^^  Moreover,  it  is  not  possible  to  bring  the 
experimental  results  of  Maccabee  and  Bell  and  the  available  theoretical 
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analyses  together  in  a  satisfactory  manner  because  of  a  number  of  conflicts 
including  the  following:  (1 )  most  of  the  results  were  in  the  form  of  peak  pressure 
measurements  which  can  be  very  misleading  because  the  peak  pressure  of  a 
wideband  signal  is  critically  dependent  on  the  response  characteristics  of  the 
receiving  system,  and  (2)  the  experimental  results  obtained  in  the  nearfield  of 
the  source  are  extremely  sensitive  to  position;  consequently,  they  were  very 
variable. 

Finally,  the  transition  from  the  linear  thermal  expansion  process  to  the 
nonlinear  evaporative  processes  appears  not  to  have  been  studied  at  all.  The 
transition  appears  to  be  caused  by  a  drastic  change  in  the  absorption  coefficient 
of  the  optical  radiation.  In  the  linear  region  of  operation,  the  absorption 
coefficient  for  1.06  pm  infrared  radiation  has  been  measured^!  to  be  in  the 
region  of  15  m-i,  but  in  the  nonlinear  region  it  is  reported^®  to  be  9  x  10^  m-T 
Without  this  change  in  the  absorption  coefficient,  the  evaporative  processes 
would  not  be  observable  at  the  currently  reported  values  of  optical  energy 
density. 

The  objective  of  this  study  is  to  bridge  some  of  these  gaps  in  the  current 
understanding  of  the  nonlinear  optoacoustic  generation  processes.  At  the 
outset  of  this  study,  it  was  proposed  that  the  work  be  approached  in  the 
following  sequence. 

(1)  The  blast  model  should  be  developed  as  far  as  possible  within  the 
linear  propagation  regime  of  the  acoustic  signal.  The  development  was  actually 
completed  under  the  previous  contract.22  The  validity  of  this  approach  will  be 
limited  to  relatively  low  laser  intensities.  The  main  advantage  in  this  approach 
is  that  it  would  allow  linear  superposition  to  be  used  in  the  reconstruction  of  the 
sound  field  and  would  avoid  problems  associated  with  nonlinear  propagation. 

(2)  The  model  will  then  be  extended  to  include  nonlinear  acoustic 
propagation  effects.  This  will  allow  the  effects  of  higher  laser  intensities  to  be 
analyzed.  The  extended  model  will  be  tested  under  controlled  conditions  to 
obtain  an  evaluation  of  its  accuracy.  The  results  obtained  will  be 


complementary  to  those  of  Maccabee  and  Bell. 

(3)  After  the  model  has  been  verified  against  experimental  results,  a 
feasibility  study  will  be  made  for  sonar  applications,  in  terms  of  energy 
conversion  efficiency,  source  levels  that  can  be  realistically  achieved,  and 
range  penetration  capabilities. 

As  it  turned  out,  the  extension  of  the  model  to  include  nonlinear 
propagation  effects  produced  a  completely  new  blast  model.  This  work  is 
described  in  Section  II.  The  model  predictions  were  compared  with 
experimental  results,  as  described  in  Section  III.  When  attempts  were  made  to 
test  the  model  at  very  high  laser  intensities,  the  experimental  results  obtained 
were  quite  unexpected.  Certainly,  they  could  not  be  explained  by  the  blast 
model.  Indeed,  the  acoustic  output  measured  was  significantly  higher  than  that 
predicted  by  the  blast  model.  Therefore,  at  this  point,  the  project  deviated  from 
the  stated  course.  A  new  optoacoustic  generation  process  was  hypothesized.  It 
was  deduced  that  an  essential  ingredient  of  the  hypothesis  is  plasma  formation. 
To  investigate  the  process  of  plasma  formation,  an  experiment  was  carried  out 
to  measure  the  photoconductivity  of  water  in  the  presence  of  t.OGp.  laser 
illumination  to  infer  the  characteristics  of  the  carriers  generated.  The  results  of 
this  experiment  are  shown  in  Section  IV.  Final  conclusions  are  drawn  in 
Section  V. 


II.  THEORY 


A  mathematical  model  of  the  explosive  reaction  that  leads  to  the 
generation  of  sound  in  water  will  be  derived.  The  reaction  is  not  self-starting.  It 
is  necessary  to  assume  that  there  exists  a  mechanism  to  absorb  the  laser 
energy  into  the  gases  that  drive  the  blast.  In  the  case  of  laser  irradiation  of  a 
solid  surface,  the  surface  must  be  sufficiently  absorptive  to  cause  local 
vaporization,  the  formation  of  a  thin  layer  of  absorptive  vapor  above  the  surface. 
This  layer  then  expands  explosively  as  it  absorbs  increasing  amounts  of  laser 
energy.  In  the  case  of  the  air-water  interface,  it  must  be  assumed  that  a 
mechanism  exists  that  allows  the  laser  energy  to  generate  a  thin  film  of  an 
absorptive  substance  from  which  the  blast  is  launched.  This  mechanism  is  not 
well  understood;  it  is  believed  to  involve  dielectric  breakdown  of  the  water 
molecules.  Let  us  for  the  moment  assume  that  a  starter  mechanism  exists  and 
proceed  to  analyze  the  basic  dynamics  of  the  blast  process.  We  will  return  to 
the  starter  mechanism  later. 

A.  Pirri's  Blast  Model 

In  the  previous  project,  a  blast  model  developed  by  Pirri'*'*  was  applied  to 
the  problem.  For  completeness,  this  model  will  be  reviewed  briefly.  The  model 
was  originally  developed  to  predict  the  impulse  imparted  by  a  laser  pulse  to  a 
solid  surface.  Pirri's  model  treats  the  explosive  reaction,  or  blast  process,  in  two 
stages  as  shown  in  Fig.  1 .  The  first  stage  is  the  one-dimensional  expansion 
stage  in  which  the  blast  expands  outward  from  the  surface.  Since  the  gases 
within  the  blast  volume  are  highly  energized,  and  since  the  blast  gases  are  in  a 
state  of  dynamic  equilibrium,  it  is  fair  to  expect  that  the  pressure  and 
temperature  within  the  expanding  blast  volume  are  approximately  uniform.  With 
this  constraint,  the  velocity  Dbo  of  the  blast  front  in  the  one-dimensional  stage 
can  be  shown  to  be  given  by^"* 

Dbo  =  [2(f-1)lo/Por  .  (1) 

and  the  pressure  Pbo  behind  the  front  is  given  by^'* 


Pbo  “  Po  ^bo  /  (  Y  +  U  > 


where  Iq  is  the  laser  intensity,  Po  is  the  density  of  the  surrounding  undisturbed 
medium  at  rest  and  y  the  ratio  of  specific  heats  of  the  blast  gas.  Using  the  gas 
equation,  the  blast  temperature  Tbo  in  °K  is  given  by 

Tb„  =  T„(  1  +  (P6<,/p„) )  ,  (3) 

where  Tg  is  the  ambient  temperature  in  °K  and  Po  the  ambient  pressure. 

The  pressure  at  the  solid  surface  Pso  is  related  to  the  pressure  behind  the 
blast  wavefront  by  a  coupling  coefficient  Cc, 

Pso  “  Pbo  , 
where  Cg  is  given  byi^ 

Cc  =  [(Y+1)/2Yf^'^^’^’  •  (5) 

The  one-dimensional  stage  terminates  when  the  height  of  the  blast  is 
approximately  equal  to  its  width.  At  this  point  the  blast  enters  the  two- 
dimensional  stage.  The  time  X20  at  which  the  changeover  takes  place  is, 
according  to  Pirri,  given  by^^ 

T2d  =  (1/D)  dh  ,  (6) 

0 

where  D  is  the  velocity  of  the  blast  front  and  d^,  the  diameter  of  the  blast  front. 
While  the  laser  beam  is  on,  D  is  equal  to  0^0,  in  which  case  X2D  'S  given  by 


Ton  =  do/^bo  . 


(7) 


Should  the  laser  pulse  cease  before  the  displacement  h  is  equal  to  do,  the  blast 
will  still  expand  in  a  one-dimensional  mode  but  at  a  declining  rate  and  the 
surface  pressure  Ps(t)  will  decline  according  to^^^ 

Ps(t)  =  Pso  (t/Xp)'^  :  for  Ta,  >t  >  Tp  .  (8) 

In  this  case,  T2d  'S  given  by 

Xjo  =  Ta,[2/(3  )1((  -  ( Dto  tp  )“1  +  tp  .  (9) 

The  second  stage  is  the  two-dimensional  expansion  phase  of  the  blast  in 
which  the  pressure  decays  rapidly  as  the  blast  area  expands.  A  cylindrically 
expanding  blast  is  assumed.  The  pressure  Ps(t)  is  given  by’'^ 

Ps{t)  =Ps(X2D)(X2D/t)  :  for  t  >  T2D  •  (10) 

As  can  be  seen  in  the  above  equations,  the  highest  pressures  are 
generated  in  the  one-dimensional  stage.  Case  studies  of  the  one-dimensional 
stage  of  Pirri's  blast  model  are  shown  in  Table  I.  Three  representative  cases 
are  shown.  In  each  case,  the  laser  pulse  is  modeled  as  having  a  peak  power  of 
4.5  MW  and  a  pulse  length  of  0.15  ps,  giving  a  total  of  energy  of  0.7  J.  The  laser 
beam  is  focused  onto  a  spot  on  a  solid  surface.  The  three  cases  shown 
correspond  to  three  different  spot  sizes,  and  hence  different  laser  intensities. 
While  the  results  shown  in  Table  I  do  not  signify  any  remarkable  differences 
between  the  three  cases,  it  will  be  shown  in  subsequent  sections  that  these 
three  cases  give  very  different  results  for  the  reaction  at  the  air-water  interface, 
and  give  some  insight  to  the  processes  involved.  In  addition,  these  cases  also 
match  the  conditions  of  the  accompanying  experimental  studies  of  laser  sound 
generation  in  water,  the  results  of  which  will  also  be  presented. 

First  of  all,  it  is  seen  that  not  all  the  laser  energy  is  used  to  drive  the  blast. 
According  to  Pirri's  theory,  once  the  blast  reaches  the  second  stage,  further 
laser  energy  arrivals  will  not  influence  the  blast  pressure  at  the  solid  surface, 
which  at  that  point  will  begin  its  decline.  Thus,  the  ratio  between  X2d  and  the 


TABLE  I 


CASE  STUDIES  OF  PIRRI'S  MODEL  FOR  THE  AIR-SOLID  INTERFACE 
CONDITIONS 


Ambient  temperature 

oc 

25 

Incident  laser  power 

W 

4.5E+6 

Pulse  duration 

s 

1.50E-07 

Energy 

J 

0.7 

CASE 

1 

2 

3 

Spot  diameter 

m 

3.2E-3 

6.6E-4 

1.7E-4 

Spot  area 

m2 

8.0E-6 

3.4E-7 

2.3E-8 

Incident  intensity 

W/m2 

5.6E+11 

1.3E+13 

2.0E+14 

Energy  density 

MJ/m2 

0.03 

1.0 

14 

AIR  BLAST  OVER  RIGID  SOLID 

Medium 

Solid 

Air 

Solid 

Air  Solid 

Air 

Molecular  weight 

gmole 

oo 

32 

OO 

32  oo 

32 

Medium  density 

kg/m3 

oo 

1.31 

oo 

1.31 

1.31 

Ratio  of  specific  heats 

N/A 

1.4 

N/A 

1 .4  N/A 

1.4 

Blast  velocity 

nVs 

N/A 

1.8E+4 

N/A 

2.6E+4  N/A 

3.6E+4 

Blast  pressure 

Pa 

5.9E+7 

1.7E+8 

1.2E+8 

3.6E-t-8  2.4E+8 

7.0E+8 

Blast  temperature 

OC 

N/A 

5.2E+5 

N/A 

1.1  E+6  N/A 

2.1  E+6 

^2D  ('*p^2Cp 

s 

N/A 

1.8E-7 

N/A 

2.6E-8  N/A 

4.7E-9 

Q 

ca 

s 

N/A 

2.1  E-7 

N/A 

2.6E-8  N/A 

4.7E-9 

Blast  front  displacement  at  T2q  m 

N/A 

3.1  E-3 

N/A 

6.6E-4  N/A 

1.7E-4 

Proportion  of  useful  energy  % 


100 


17,2 


3.16 


laser  pulse  length  Xp  may  be  taken  as  an  approximate  indication  of  the  useful 
fraction  of  the  laser  energy.  It  is  seen  that,  while  100%  of  the  energy  is  utilized 
in  case  1,  in  cases  2  and  3  only  small  percentages  of  the  laser  energy,  17.2% 
and  3.16%,  respectively,  are  effectively  utilized.  This  is  because,  in  cases  2  and 
3,  the  duration  of  the  first  stage  of  the  blast  X2d  is  very  short  compared  to  the 
laser  pulse  length.  The  computed  blast  pressures,  velocities,  and  temperatures 
are  also  shown.  It  can  be  seen  that  the  computed  blast  velocities  are  all 
supersonic.  The  blast  pressures  and  the  pressures  transmitted  to  the  solid 
surface  are  all  estimated  to  be  in  excess  of  100  atmospheres  (1  atmosphere  = 
10®  Pa).  The  values  are  most  extreme  in  the  third  case,  where  the  laser  spoL  is 
smallest  and  the  laser  intensity  is  highest.  The  computed  temperatures  of  the 
blast  gases  range  from  0.52  to  2.1  million  °C. 

In  a  preliminary  study,23  this  model  was  applied  with  little  modification  to 
the  optoacoustic  sound  generation  process.  It  was  assumed  that,  for  the  most 
significant  part  of  the  blast,  the  water  surface  remained  unperturbed  and 
therefore,  for  practical  purposes,  behaved  as  a  solid  surface.  The  surface 
pressure  that  the  model  predicted  would  form  a  distributed  source  function  of 
the  underwater  acoustic  signal.  The  acoustic  signal  was  then  constructed  by 
linear  superposition  of  the  Green's  function  from  each  elemental  source  point. 
Although  it  was  a  very  crude  model,  it  gave  very  reasonable  results.  The 
following  is  a  more  rigorous  treatment  of  the  process  of  optoacoustic  sound 
generation,  starting  from  the  basic  hydrodynamic  principles. 

B.  Blast  Model  Development  for  the  Air-Water  Interface 

To  overcome  the  shortcomings  of  the  Pirri  model,  a  new  mode!  will  be 
derived  from  the  basic  hydrodynamic  equations.  Let  us  start  by  considering  the 
very  beginning  of  an  optically  induced  reaction  at  the  air-water  interface.  At  the 
very  start  of  the  reaction,  any  activity  would  be  confined  to  a  thin  layer  whose 
vertical  dimension  is  expected  to  be  on  the  order  of  a  wavelength  of  the  laser 
light  or  the  size  of  a  water  molecule,  whichever  is  greater.  In  comparison,  the 
laser  beam  would  be  much  wider;  therefore,  the  initial  process  may  be 


considered  as  a  one-dimensional  process.  It  is  postulated  that  there  is  a 
pressure  buildup  due  to  thermal  absorption.  Therefore,  there  must  be  a 
pressure  discontinuity  propagating  into  the  air  counter-balanced  by  a  pressure 
front  propagating  into  the  water.  Initially,  this  may  be  considered  as  a  one¬ 
dimensional  process.  When  the  vertical  separation  of  the  two  fronts  reaches  a 
value  where  it  begins  to  significantly  affect  the  course  of  the  reaction,  then  the 
one-dimensional  assumption  will  break  down.  The  division  of  the  evolution  into 
one-  and  two-dimensional  stages  is  a  useful  simplification  that  is  helpful  in 
gaining  an  insight  into  the  problem. 

The  new  model  for  the  one-dimensional  stage,  as  shown  in  Fig.  2,  shows  a 
discontinuity  front  traveling  into  the  air  and  another  into  the  water  from  the 
interface.  For  equilibrium,  the  pressures  behind  the  two  fronts  must  be  equal. 
At  this  point,  the  two  separate  fronts  are  launched  outward  from  the  interface. 
The  fronts  will  be  called  blast  fronts  or  waves,  since  they  are  driven  by  hot 
gases  which  are  energized  by  the  laser  beam. 

First  of  all,  consider  the  equations  governing  each  front.  The 
hydrodynamic  equations  governing  a  well  defined  planar  discontinuity  can  be 
found  in  a  number  of  well-known  texts.^^*^^  They  are  as  follows. 

The  conservation  of  mass  equation, 

pu  =  PoD,  (11) 

the  momentum  equation. 


p  +  pu2  =  PqD^, 

and  the  energy  equation, 

lo/(PoD)  =  CO  +  0.5u2  -  0.5D2  , 


(12) 


(13) 
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where  p  is  the  excess  pressure  above  the  ambient  pressure  Po,  p  and  u  are  the 
density  and  particle  velocity  of  the  gas  behind  the  front,  u  is  measured  relative 
to  the  moving  front,  Iq  is  the  power  supplied  to  the  planar  discontinuity  per  unit 
area,  po  is  the  density  of  the  unperturbed  fluid  (gas  or  liquid)  which  is  assumed 
to  be  at  rest,  D  is  the  velocity  of  the  front  relative  to  the  unperturbed  fluid,  and  co 
is  the  specific  enthalpy  increase. 

1.  Maximum  Entropy  Condition 

For  a  gas,  the  specific  enthalpy  to  increase  is  of  the  form 

0)  =  (R/m)(T-To)7/((Y-1) 

=  (P  +  Po)Y/(  (Y  - 1  )p)  *  PoY/(  (Y  - 1  )Po).  (1 4) 

where  T  is  the  gas  temperature,  Tq  and  Po  the  ambient  temperature  and 
pressure,  R  the  universal  gas  constant,  and  m  the  molecular  weight. 

If  it  can  be  assumed  that  p»Po,  then  it  is  seen  that  Po  may  be  neglected. 
Eliminating  D  and  u  from  the  above  equations,  Raizer^^  showed  that  a  simple 
solution  of  the  pressure  may  be  found  in  terms  of  the  laser  intensity  Iq.  the 
physical  constants  po  and  y.  and  the  expansion  coefficient  rj.  Giving  the 
pressure  term  p  subscripts  "bo"  to  indicate  that  it  pertains  to  the  pressure  behind 
the  lilast  front  in  the  Qne-dimensional  stage,  the  solution  is  given  by 

Pbo  -  Po’'^ io)/(ti(ni)/(r-i)  -  1 )  .  (15) 


where  ti  is  defined  as 


The  expansion  coefficient  ti  needs  to  be  addressed  at  some  length.  There 
is  nothing  in  the  above  equations  that  can  determine  the  value  of  ti.  According 
to  Raizer23  and  Cole,^"^  the  value  is  determined  by  the  Chapman-Jouguet 
condition.  While  Chapman25  simply  assumed  that  the  front  must  travel  at  the 
minimum  velocity  consistent  with  the  above  equations,  Jouguet^^  postulated 
that  the  front  must  travel  at  a  velocity  equal  to  Cq  +  (D-u),  that  is,  the  velocity  of 
small  amplitude  waves  Cq  plus  the  absolute  particle  velocity  (D-u).  Both 
conditions  are  equivalent  and  lead  to  identical  results.  There  are  a  number  of 
physical  arguments  that  can  be  made  in  support  of  the  Chapman-Jouguet 
condition. 

In  general,  the  value  of  ti  may  be  explained  by  the  principle  of  increasing 
entropy,  a  corollary  of  the  second  law  of  thermodynamics.  Since  the  front  is  a 
discontinuity,  the  thermodynamic  state  of  the  fluid,  as  it  passes  through  the 
discontinuity,  must  for  a  very  brief  moment  be  indeterminate.  That  is,  it  can  go 
through  the  full  range  of  possible  states.  Since  the  principle  of  increasing 
entropy  only  allows  transitions  to  states  with  higher  entropy,  it  follows  that  the 
final  state  will  be  the  one  with  the  maximum  entropy.  Therefore,  the  fluid 
emerges  from  the  discontinuity  at  the  maximum  entropy  allowed  by  the 
governing  equations. 

The  maximum  entropy  condition  may  be  simply  found  by  taking  the  partial 
derivative  of  the  entropy  S  with  respect  to  tj  and  setting  dS/dr]  equal  to  zero. 
For  a  gas  initially  at  pressure  Po  and  density  Po,  and  arriving  at  the  state  given 
by  P+Po  and  t|,  the  entropy  increase  may  be  computed  by  considering  the 
transition  to  be  in  two  stages:  an  isentropic  expansion  to  Po'  and  r\  followed  by  a 
constant  volume  heating  to  the  final  state.  From  the  adiabatic  gas  equation,  the 
pressure  po'  is  given  by 

Po’  =  Pon"^  (17) 


at  an  absolute  temperature  Tq  of 


To'  =  ( Po'Tl/Po  )(m/R) 


(18) 


where  m  is  the  molecular  weight  of  the  gas.  and  R  the  universal  gas  constant.  In 
heating  at  constant  volume  from  pressure  Po*  to  p+Po,  the  final  temperature  T  is 
given  by 

T  =  ( (Pbo+Po)Tl/Po  )(m/R)  .  (19) 

The  entropy  increase  S  is  given  by 

S  =  c,/TdT 

’o 

=  c,  ( ln(T)  -  In(To-) )  .  (20) 

where  Cy  is  the  specific  heat  at  constant  volume.  Substituting  for  Tg'  and  T  from 
the  above  equations  and  differentiating  with  respect  to  -n,  the  result  is 

as/aii  =  Cv(  (ai/an)  /t)  -  (aiovari)  rr^') ) 

=  (Cv/h)(Y  +  (3p/^)Tl/(Pbo+Po) )  •  (21) 

Substituting  for  p^o  from  Eq.  (15),  carrying  out  the  partial  derivative  ap/ari,  and 
then  setting  as/a'n=0,  the  value  of  t]  at  the  maximum  entropy  condition  may  be 
found.  Cole^"^  has  shown  that,  where  the  pressure  is  large  compared  to  the 
atmospheric  pressure,  the  condition  as/aq  =  0  is  equivalent  to  the  Chapman- 
Jouguet  condition.  Thus,  the  Chapman-Jouguet  condition  may  be  justified  on 
the  basis  of  the  laws  of  thermodynamics.  According  to  Raizer,  the  Chapman- 
Jouguet  condition  for  a  gas  occurs  when 


n  =  Y/(Y+1)  • 


(22) 


Let  us  now  return  to  the  fundamental  equations,  Eqs.  (11)-(13),  and 
derive  the  solutions  for  the  air-water  interface  without  neglecting  the  ambient 
pressure  term.  Water,  however,  is  far  from  being  a  perfect  gas;  therefore,  it  is 
necessary  to  deal  with  certain  complications  in  the  formulation  of  its  pressure- 
density-temperature  relationships. 

It  is  not  possible  to  give  a  complete  description  of  the  thermodynamic 
behavior  of  water  through  simple  analytical  formulae.  For  the  purposes  of  this 
study,  we  will  consider  two  idealized  regimes,  the  gas  regime  and  the  liquid 
regime,  which  will  be  defined  below.  These  two  regimes  were  found  to  be 
adequate  most  of  the  time.  Beyond  and  between  these  two  regimes,  there  is  a 
vast  area  of  state  space  that  will  be  dealt  with  only  as  the  need  arises.  The 
criteria  for  the  gas  regime  were  found  empirically  by  searching  the  state  space 
as  contained  in  the  steam  tables.^®  It  was  found  that  in  certain  regions  of  the 
state  space,  steam  behaved  very  much  like  a  perfect  gas,  i.e.,  its  properties 
differed  little  from  that  of  a  perfect  gas  of  the  same  molecular  weight.  The 
idealized  liquid  regime  is  an  adaption  of  the  state  equations  for  high 
temperatures  given  by  Cole.  The  criteria  for  the  liquid  regime  were  also  found 
empirically  by  searching  the  steam  tables.  The  pressure  and  temperature 
ranges  of  the  available  steam  tables  are  somewhat  limited,  and  beyond  these 
ranges  one  can  only  extrapolate. 

In  the  following,  the  equations  for  the  pressure  and  front  velocity  will 
be  rederived  without  neglecting  the  ambient  pressure  terms,  and  the  analysis  is 
extended  to  include  the  possibility  of  superheated  steam  or  liquid  water  as  a 
blast  fluid. 


a.  Gas  Regime 


In  the  gas  regime,  all  of  the  water  absorbed  by  the  blast  is 
converted  into  superheated  steam.  The  blast  absorbs  water  at  the  rate  given  by 


the  velocity  of  the  underwater  front.  In  this  regime,  the  specific  enthalpy  term  in 
Eq.  (13)  may  be  approximated  by 


CO  —  cog  —  (R/rn3)(Tg  -  Tg(^)Ys/(Ys ' 

=  (Pso +Po)Ys/((Ys- 1)Ps)  *  PoYs/({Ys- 1)Pow)  +  Esh  ■  (23) 

All  the  terms  that  pertain  to  steam  have  been  given  a  subscript  "s”  and  those 
that  pertain  to  water  the  subscript  "w",  Pgo  is  the  steam  pressure  behind  the  front 
in  the  one-dimensional  stage,  Tg  is  the  temperature,  pow  is  fhe  density  of  the 
unperturbed  water,  and  mg  is  the  molecular  weight  of  steam,  which  is  18.  The 
term  egh  is  the  latent  of  vaporization  at  ambient  density  pow  and  pressure  Po,  and 
Tgh  is  the  interpolated  temperature  at  ambient  pressure  and  density  of  the 
undisturbed  water,  as  predicted  by  the  gas  equation.  The  term  Egh  is  given,  in 
terms  of  the  latent  heat  at  constant  volume  and  temperature  Eg^,  by 


^  ~  Sgy  -  Po(PoW^Psv)  YsA  (Ys  *  ^  )Pow) 

=  Esv  -  3300  kJ/kg  ;  at  25'’C  and  atmospheric  pressure  ,  (24) 

where  psv  is  the  interpolated  density  of  steam  at  ambient  temperature  Tq  and 
pressure  Po  as  determined  by  the  gas  equation. 

These  equations  may  also  be  applied  to  a  perfect  gas  by  simply  setting  the 
latent  heat  term  Egh  equal  to  zero,  and  Tgh  equal  to  T^.  Thus  it  may  also  be 
applied  to  the  air  blast  propagating  upwards  from  the  interface.  To  distinguish 
the  two  cases,  the  subscript  "a"  will  replace  "s"  or  "w"  for  the  terms  associated 
with  air  blast.  For  air  the  effective  molecular  weight  mg  is  29. 

The  ratio  of  specific  heats  y  for  many  gases,  including  air,  is  dependent  on 
temperature  and  very  weakly  dependent  on  pressures.  For  air,^^  the  value  of  Ya 
varies  between  1 .4  at  room  temperature  to  1.23  at  2800°K.  For  steam,  however, 
the  variations  are  much  larger.  Using  steam  tables,^®  the  values  of  Ys  were 


computed  over  the  available  range  of  temperatures  and  pressures  as  shown  in 
Fig.  3.  It  can  be  seen  that  for  temperatures  above  600°C  the  value  of  Ys  lies 
between  1 .2  and  1 .3.  Therefore,  a  compromise  value  of  1 .25  was  chosen  for  Ys- 
This  is  comparable  to  the  y  values  of  other  triatomic  gases  such  as  SO2  at  1 .26 
and  CO2  at  1.31. 

The  latent  heat  of  vaporization  at  a  constant  volume  is  the  difference 
between  the  internal  energies  of  water  and  steam.  In  Fig.  4,  plots  of  the  internal 
energies  of  water  and  as  a  function  of  temperature  at  various  pressures, 
obtained  from  the  steam  tables,  are  compared  with  a  straight  line  perfect  gas 
law  computed  for  the  molecular  weight  of  water  and  a  specific  heat  ratio  of  1 .25. 
It  is  seen  that  at  high  temperatures,  the  curves  tend  towards  the  perfect  gas  line. 
Assuming  an  initial  state  at  25°C  and  atmospheric  pressure,  it  is  seen  that  the 
latent  heat  is  approximately  2300  kj/kg.  Substituting  into  Eq.  (24),  this  gives 
Esh  a  value  of  -1000  kJ/kg. 


In  addition  to  the  gas  regime,  the  data  in  Fig.  4  also  suggest  a 
liquid  regime.  In  the  liquid  regime,  the  specific  enthalpy  of  water,  adapted  from 
the  state  equations  given  by  Cole,^^  is  defined  by  the  relationship 

0)  =  =  ((Pwo  Po)Yw  (Pow  ■  Pw)V((Yw  ■  ^  )pw)  ■  PoYv/({Yw  ■  ^  )Pow) 

=  Cwv(Tw  -  Tq)  +  (Pyyo  +  Po)/pw  ■  Pc/Pow  •  (25) 

where  all  the  terms  that  pertain  to  the  liquid  regime  have  been  given  the 
subscript  "w",  p^^o  is  the  excess  pressure  of  liquid  water  above  the  ambient 
pressure  po  in  the  one-dimensional  stage,  c^^  is  the  specific  heat  at  constant 
volume  which  is  approximately  a  constant,  T,^  is  the  absolute  temperature 
behind  the  front  and  Tq  is  the  ambient  temperature  in  front  of  it,  which  is 
assumed  to  be  298°K  (25°C)  for  calculation  purposes,  and  c^o  is  the  speed  of 
small  amplitude  sound  waves  in  the  unperturbed  water  which  at  25°C  is  equal 
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FIGURE  3 

RATIO  OF  SPECIFIC  HEATS  OF  STEAM  COMPUTED  FROM  STEAM  TABLES 
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to  1534.6  m/s.  Using  values  from  the  steam  tables,  it  was  estimated  that  a 
reasonable  approximation  may  be  obtained  with  Yw  =  l  Pow  =  ‘•001  kg/m^, 
and  Cwv  =  4.2  kJ/kg.  The  adiabatic  state  equation  is  given  by^^ 


Pwo  -  Pow/Tw)(  (Pw^Pow)^'*^  *  ^)  ■  (2®) 

It  is  interesting  to  note  that  if  the  asymptotic  functions  of  the  internal  energy 
of  the  liquid  and  gas  regimes,  as  shown  in  Fig.  4,  were  extrapolated,  they  would 
intersect.  There  is  a  corresponding  intersection  of  the  enthalpy  functions,  which 
is  found  by  equating  the  enthalpies  as  given  by  Eqs.  (23)  and  (25);  the 
temperature  at  the  intersection  is  found  to  be  978°C.  Beyond  this  temperature, 
the  liquid  regime,  if  it  existed,  would  apparently  have  a  higher  specific  enthalpy 
than  the  gas  regime. 

Despite  appearances,  Eq.  (25)  is  not  a  gas  equation,  and  Yw  not  the  ratio 
of  specific  heats.  The  specific  heat  at  constant  volume  Cwv,  as  calculated  from 
the  steam  tables,  was  found  to  be  approximately  constant  in  this  regime. 

The  regions  on  the  pressure-temperature  diagram  in  which  the  above 
equations  for  the  gas  and  liquid  regimes  are  valid  are  shown  in  Fig.  5;  in  these 
regions  the  parameter  estimates  of  the  two  respective  models  were  found  to  be 
within  1 0%  of  the  tabulated  values  from  steam  tables.  Beyond  the  range  of 
tabulated  values,  the  validity  of  the  models  is  not  known  with  any  certainty:  the 
boundaries  shown  are  mere  extrapolations  of  the  boundary  functions  found  for 
the  tabulated  region. 

c.  Blast  Pressure  Solution 

Eliminating  D  and  u  from  Eqs.  (11)-(13),  a  general  relationship 
between  enthalpy,  pressure,  and  coefficient  of  expansion  is  obtained. 


PRESSURE  -  Pa 


LIQUID 

REGIME 


GAS 

REGIME 


LIQUID 

REGIME 


AN[i 

SATURATED 
STE/iM  i 


r  GAS 
REGIME 


REGIIOIM 


TEMPERATURE  -  C 


FRONT 

STATES 


INTERIOR 

STATES 


OCase  1 
O  Case  2 


^Case  1 
^  Case  2 


I  Case  3 


,  Case  3 


FIGURE  5 

LIQUID  AND  GAS  REGIMES  OF  THE  STATES  OF  WATER 


ARL;UT 
AS-88-83 
NPC  -  GA 
1-27-88 


Substituting  for  co  from  either  Eq.  (23)  or  (25).  and  applying  subscripts  "w"  or  "s" 
where  appropriate,  a  cubic  equation  in  p^^^  is  obtained,  of  the  form 


Pxo^  +  SaxPxo^'^  +  2bx  =  0 


For  the  pressure  behind  the  one-dimensional  front  in  air  Pao,  the  coefficients  ag 
and  ba  are  given  by 

a,  =  (2/3)(  Po7a(na  -  1  )/(Ya  '  '  M  W  Ya  +  1  )/(Ya  -  1 ))  '  ' )  . 

ba  =  -loa(1-na)'®Poa''^/((nd(Ya+iy(Ya  '  1))  '  1)  ■ 

where  log  is  the  power  density  absorbed  by  the  air  front. 

Assuming  that  all  the  absorbed  water  is  changed  into  steam,  the  pressure 
behind  the  underwater  front  Puo  'S  given  by  the  steam  pressure  Pso-  The 
coefficients  and  bg  for  solving  for  Pgo  are  given  by 

ag  =  (2/3)[  Egh  pow  +  YsPo{Tls  - 1  )/(Ys  - 1  )H(ns(  Ys  +  1  )/(Ys  *  1 ))  *  1 1 . 

bs  =  *lou  Pow^^/((^s(Ys  +  1V(Ys  -  1))  -  1)  . 

where  Iqu  is  the  power  density  at  the  underwater  front. 

Finally,  assuming  that  all  the  absorbed  water  remains  in  the  liquid  state, 
the  pressure  behind  the  underwater  front  p^o  is  given  by  the  water  pressure 
PwQ.  The  coefficients  a^  and  by^  for  solving  for  Py^^  are  given  by 

a^  =  (2/3)po^(c^  +  YwPc/Pow)(  ^ V{  ^w(  Yw  ^ )  *  ( Yw "  > 

bw  “  '^OU  (^*^w)  Pow  A  (T1w(  Yw  ^  )/(Yw  ■  f))  ■  f)  • 


The  general  solution  where  a  valid  one  exists  is  given  by 

Pxo  =  [  (-bx  +  ''(bx"  +  ax“) )“  -  (bx  +  '/(bTT^) )’«  f  .  (29) 

where  the  condition  for  a  valid  solution  is:  b)(^+ax^>0.  This  condition  implies 
that  there  are  minimum  levels  of  the  power  densities  or  intensities  log  and  Iqu 
that  are  necessary  to  the  sustenance  of  the  respective  fronts. 

In  each  case,  the  value  of  h^  remains  to  be  found  by  the  maximum  entropy 
condition.  For  a  gas  at  high  pressure,  including  air  and  superheated  steam,  it  is 
given  by  the  Chapman-Jouguet  condition  in  Eq.  (22).  In  the  general  case, 
however,  it  must  be  found  by  the  maximum  entropy  condition.  Since  in  all  three 
cases  that  are  of  interest  here  the  specific  heats  at  constant  volume  may  be 
treated  as  constants,  the  procedure  for  obtaining  an  expression  for  the  entropy 
in  terms  of  r[^  is  straightforward.  Using  the  state  equation  such  as  Eqs.  (14), 
(23),  or  (25)  as  the  case  may  be,  the  fluid  is  allowed  to  expand  isentropically  to 
an  intermediate  state  that  has  the  same  density  as  its  final  density.  An 
expression  is  found  for  the  temperature  at  this  intermediate  state.  Then,  it  is 
allowed  to  heat  up  at  constant  volume  to  its  final  state,  and  a  second  expression 
is  found  for  the  temperature  at  the  final  state.  Since  the  specific  heat  at  constant 
volume  is  approximately  constant,  the  result  is  an  expression  similar  to  Eq.  (20). 
Finally,  the  expression  for  entropy  is  differentiated  with  respect  to  to  find  the 
maximum  entropy  condition. 

For  air  and  superheated  steam  at  high  pressures,  the  equations  for  the 
variable  part  of  the  entropy  are  identical  to  those  of  the  air  blast,  and  thus 
Eqs.  (17)-(21)  may  be  used. 

For  the  liquid  regime,  the  expressions  for  the  intermediate  and  final 
absolute  temperatures  Ty^,o'  and  T^^,  corresponding  to  T^'  and  T,  are  not  as 
simple  as  for  the  gas  regime.  They  are  found  from  the  state  equation  of  Eq.  (25) 
and  the  adiabatic  equation  of  Eq.  (26)  to  be 


“  ((^v/9w)(^w  V(9w’^  )  *  Po(^  "^wV{^ow(9w''^ )  )  '  (30) 

Cwv  (Tyy  -  Tq) 

~  (Pwo^w  ■*■  *'ow(^w‘^  ))A*’ow(9w‘^ ))  *  Po(^  ■^w)^(*^Ow(9w'^  )  )  (31) 

Substituting  for  Tq'  and  T  in  Eq.  (20)  with  the  above  expressions  for  Ty^o'  and  T^, 
and  differentiating  with  respect  to  hy,,,  the  maximum  entropy  condition  may  be 
found.  The  resulting  equation,  however,  is  very  messy.  Thus,  for  this  case,  a 
numerical  method  for  finding  the  maximum  entropy  condition  may  be 
preferable. 

Numerical  results  for  the  three  case  studies  are  shown  in  Table  II.  The 
maximum  entropy  conditions  for  all  the  blast  fluids  were  found  numerically  by 
the  Newton  interpolation  method.  The  values  of  h  for  maximum  entropy  and  for 
the  Jouguet  condition  are  shown  for  comparison.  It  is  seen  that,  where 
applicable,  they  are  identical,  except  in  case  1  where  the  pressure,  at 
130  atmospheres  (1.3  ¥  10^  Pa),  is  sufficiently  low  that  atmospheric  pressure  is 
not  entirely  negligible.  The  minimum  intensity  required  to  sustain  each  front, 
calculated  according  to  the  validity  condition  for  Eq.  (29),  is  shown  along  with 
the  actual  power  density  absorbed;  in  case  1,  it  will  be  shown  that  the  laser 
intensity  is  insufficient  to  support  the  total  process,  and  therefore  the  minimum 
intensity  is  used  for  calculation  purposes.  In  all  cases,  it  was  found  that  the 
water  exiting  the  underwater  front  is  in  the  liquid  regime,  and  the  air  at  the  air 
front  remains  in  the  gaseous  phase.  It  should  be  noted  that  the  calculated  air 
temperature  as  it  exits  the  air  front  is  very  high,  ranging  from  22,000°C  in 
case  1  to  1.9  million°C  in  case  3.  It  will  be  shown  in  the  following  section  that 
the  air  is  not  expected  to  remain  at  these  extremely  high  temperatures  for  any 
significant  length  of  time  because  it  will  immediately  mix  with  the  gases  in  the 
interior  of  the  blast  which  are  considerably  cooler. 


TABLE  II 

EQUILIBRIUM  BLAST  MODEL  CASE  STUDIES.  FRONT  CONDITIONS 


PHYSICAL  CONSTANTS 


Medium 

Steam 

Water 

Air 

Molecular  weight 

g  mole 

18 

18 

32 

Ambient  density 

kg/m3 

N/A 

1001 

1.3 

Adiabatic  exponent 

1.25 

1.335 

1.4 

Latent  heat  of  vaporizatton 

J/kg 

N/A 

2.3E+6 

N/A 

Specific  heat  at  constant  volume  kJ/kg 

4.2 

4.2 

0.72 

Sound  speed 

m/s 

N/A 

1535 

330 

CASE 

1 

2 

3 

Spot  diameter 

m 

0.0032 

0.00066 

0.00017 

Spot  area 

m2 

8.0E-6 

3.4E-7 

2.3E-8 

Incident  intensity 

W/m2 

5.6E+11 

1.3E+13 

2.0E+14 

Energy  density 

MJ/m2 

0.03 

1.0 

14 

FRONT  CHARACTERISTICS 

Medium 

Water 

Air 

Water 

Air 

Water 

Air 

Minimum  intensity 

W/m2  1.7E+10 

2.6E+7 

1.7E+10 

2.6E+7 

1.7E+10 

2.6E+7 

Actual  intensity 

W/m2  1.7E+10 

7.7E+10 

1.0E+12 

2.5E+12 

2.2E+13 

6.3E+13 

Ti  ( maxirrxim  entropy) 

0.9960 

0.5852 

0.9856 

0.5833 

0.9780 

0.5833 

11  (Jouguet  condition) 

N/A 

0.5833 

N/A 

0.5833 

N/A 

0.5833 

Phase 

LIQUID 

GAS 

LIQUID 

GAS 

LIQUID 

GAS 

Front  velocity 

m/s  1 .8E+3 

4.9E+3 

3.0E+3 

1.5E+4 

7.1  E+3 

4.5E+4 

Mach  number 

1.16 

14.71 

1.96 

46.81 

4.64 

137.13 

Absolute  particle  velocity 

m/s  7.2 

2.0E+3 

43 

6.4E+3 

1.6E+2 

1.9E+4 

Front  pressure 

Pa  1.3E+7 

1.3E+7 

1.3E+8 

1 .3E+8 

1.1  E+9 

1.1  E+9 

Front  temperature 

OC  27 

2.2E+4 

92 

2.3E+5 

7.6E+2 

1.9E+6 

Specific  enthalpy 

kJ/kg  2.2E+4 

2.0E+7 

4.1E+5 

2.0E+8 

4.2E+6 

1.7E+9 

Specific  entropy  kJ/kg/K  32 

2.9E+3 

8.5E+2 

4.6E+3 

5.2E+3 

6.1  E+3 

Notes  : 

1 .  Ambient  temperature,  laser  power,  and  pulse  duration  are  as  stated  in  Table  I. 

2.  The  intensity  at  the  underwater  tront  in  case  1  (in  italics)  is  the  minimum  consistent  with  the 
maximum  entropy  condition. 


To  obtain  equilibrium,  the  pressures  behind  the  air  and  the 
underwater  fronts  Pao  and  p^o.  as  defined  by  Eq.  (28)  and  its  solution  Eq.  (29), 
and  the  pressure  Pgo  of  the  interior  of  the  blast  volume  must  all  be  equal.  Thus, 

Pgo  =  Pao  =  Puo  •  (32) 

With  two  fronts  back-to-back  producing  gases  at  the  same  pressure  but  at 
different  temperatures  and  densities,  there  is  bound  to  be  an  increase  in 
entropy  in  the  subsequent  temperature  equalization.  This  will  result  in  a  net 
loss  of  pressure  in  the  interior  of  the  blast.  To  maintain  equilibrium,  a  proportion 
of  the  supplied  laser  power  must  be  diverted  to  the  gases  in  the  interior  in  order 
to  restore  the  gas  pressure. 

The  initial  enthalpy  density  rate  hg'  entering  the  interior  per  unit  area  of  the 
fronts  is  given  by 

hg  *=  (^ePaoPoa  ^u^uoPou)  >  (33) 

where  Dgo,  Duo.  Poa-  and  Pou  are  the  front  velocities  in  air  and  underwater  and 
the  unperturbed  densities  of  the  two  media,  respectively.  The  specific 
enthalpies  (Og  and  (0^  may  be  computed  from  Eqs.  (23)  and  (25),  respectively. 
The  final  enthalpy  density  rate,  after  temperature  equalization,  may  be 
expressed  in  a  similar  form, 

hg  =  (<j^ag^aoPoa  *^ug^uoPou)  >  (34) 

where  cOgg  and  co^g  are  the  final  specific  enthalpies  after  temperature 
equalization  of  the  fluids  in  the  interior.  The  subscript  "g"  is  added  to  distinguish 
the  physical  properties  at  equilibrium  from  the  corresponding  properties  of  the 
fluids  in  their  initial  conditions.  The  energy  absorption  rate  log  of  the  blast 
interior  per  unit  frontal  area  is  then  given  by 


The  terms  cogand  0)^  may  be  evaluated  through  Eq.  (23).  To  find  cogg  and 
o)ug,  we  start  by  finding  the  corresponding  final  component  densities  pug  and  pgg 
of  the  interior  fluids.  One  equation  relating  them  is  obtained  by  invoking 
conservation  of  mass  and  total  volume,  giving 

(Dgo  +  Duo)  =  (DaoPoa^Pag)  (DyoPou^Pug)  •  (^®) 

A  second  independent  expression  relating  pgg  and  pug  nnay  be  found  by 
equating  the  temperatures  of  the  two  fluids  and  setting  their  pressures  equal  to 
the  blast  pressure  Pg,  through  the  state  equations  of  Eqs.  (23)  and  (25).  Thus, 
it  is  possible  to  solve  for  pgg  and  pug- 

For  the  case  where  all  the  absorbed  water  is  converted  into  superheated 
steam,  it  is  found  that  the  density  ratios  at  equilibrium  are  simply  given  by  the 
molecular  weight  ratios,  giving 

Pag/pug  =  Hl^mg  .  (37) 

Using  Eqs.  (36)  and  (37),  the  component  densities  pgg  and  pug  may  be  simply 
solved  giving,  for  pug. 

Pug  “  ((DaoPoa  DuoPou)ADao  ■*"  D^q)  (38) 

and  a  symmetrical  expression  for  pgg.  Setting  the  pressures  equal  to  Pg^  and 
using  the  above  expression  for  pgg  and  p^g,  expressions  for  the  final  specific 
enthalpies  cogg  and  coyg  of  the  two  gaseous  components  may  be  obtained 
through  the  state  equation  of  Eq.  (23).  Finally,  substituting  the  resulting 
expressions  for  cOgg  and  (o^g  into  Eq.  (34),  a  solution  for  the  final  enthalpy 
density  rate  may  be  found  for  the  air-steam  mixture. 


In  the  case  where  the  absorbed  water  remains  a  liquid,  the  problem  is 
more  complicated.  There  are  three  distinct  possibilities  here.  One  is  that  the 
liquid  water  within  the  blast  may  coalesce  into  large  droplets  resulting  in  very 
little  heat  exchange  with  the  surrounding  hot  air,  in  which  case  there  will  be  no 
temperature  equalization.  This  possibility  is  considered  unlikely  since  the 
gases  from  the  air  front  will  be  at  a  high  temperature  and  highly  agitated  and 
therefore  it  would  tend  to  break  up  the  water  droplets.  The  second  possibility  is 
that  the  liquid  water  evaporates  into  steam  on  contact  with  the  hot  air;  in  this 
case  Eqs.  (33)-(38)  would  be  applicable.  The  third  is  that  the  liquid  water  exists 
as  a  fine  spray  and  exchanges  heat  with  the  hot  air  but  remains  in  liquid  form. 
This  last  possibility  is  addressed  by  equating  the  temperatures  of  the  air  and  the 
liquid  water  and  setting  the  pressures  equal  to  equilibrium  pressure  Pg. 
Through  the  state  equations  of  Eq.  (23)  and  (25),  the  second  independent 
equation  relating  Pag  and  pug  is  given  as 

((Pg  Po)yw  ^wo  (Pow  *  Pug)  )A(Yw  "  ^  )Pug)  *  PoYw/({Yw  ‘  ^  )Pow) 

=  Cwv(  (Pgl^a  A^Pag))  '  "^o)  "*■  (Pg  +  Po)^Pug  *  Po^Pow  •  (39) 

Eliminating  Pag  from  Eqs.  (36)  and  (39),  pgg  is  given  by 

Pug^Pow  ~  (  (Yw  *  ^  )^^Pg(^a^^)(^uo*'ou^^aoPoa)  Pg  ■*"  Pow  Po) 

((Yw  '  "I  )(^^vvPg(nia /R)(Duo  +  Dao)/(DaoPoa)  '  "^o)  Pow  ■*"  Po)  •  (^9) 


Substituting  for  r^g  from  Eq.  (40)  into  (36)  a  similar  expression  for  pag  may  be 
found.  Applying  the  expressions  for  the  final  densities  of  the  air  and  water 
components,  pag  and  p^g,  to  their  respective  state  equations,  Eqs.  (23)  and  (25) 
respectively,  expressions  for  the  final  specific  enthalpies  of  the  two  components 
may  be  obtained.  Again,  substituting  the  resulting  expressions  for  coag  and  cOyg 
into  Eq.  (34),  a  solution  for  the  final  enthalpy  density  rate  may  be  found  for  the 
air-water  mixture. 


From  the  above  analysis,  it  is  seen  that  the  laser  intensity  Itot  niust  be 
divided  into  three  parts, 

^ot  “  ^oa  ^ou  ^og  • 

where  loa  is  the  part  supplied  to  the  air  front,  Iqu  to  the  underwater  front,  and  log 
to  the  interior  of  the  blast.  The  pressure  equilibrium  conditions  as  stated  in 
Eq.  (32)  determine  the  exact  division  of  intensities. 

Numerical  case  studies  of  the  equilibrium  conditions  are  shown  in 
Table  III.  The  values  of  log,  Iqu-  log  were  found  iteratively,  using  the 
pressure  and  temperature  equilibrium  conditions  as  defined  in  the  above 
equations.  The  percentage  breakdown  of  the  power  distribution  shows  that  a 
very  large  proportion  is  absorbed  by  the  interior  of  the  blast.  In  case  1 ,  we  have 
a  situation  where  the  laser  intensity  is  sufficient  to  drive  the  two  fronts,  but 
insufficient  to  maintain  the  internal  pressure  of  the  blast;  therefore,  the  blast  is 
effectively  quenched  as  soon  as  the  fluids  from  the  two  fronts  meet.  It  is  also 
apparent  that  only  a  proportion  of  the  laser  intensity,  6.39%  in  case  2  and 
11.10%  in  case  3,  is  directly  employed  in  driving  the  underwater  blast  front. 
These  percentages  are  expected  to  have  a  direct  influence  on  the  sound 
generation  efficiency  of  the  process.  In  cases  2  and  3,  it  was  found  that  the 
water,  which  was  still  in  a  liquid  form  as  it  emerged  from  the  underwater  front, 
would  be  converted  into  superheated  steam  (gas  regime)  after  mixing  with  the 
hot  air  from  the  air  front;  the  temperature  of  the  mixture,  while  still  high,  at 
1,747°C  in  case  2  and  17,940°C  in  case  3,  is  still  significantly  lower  than  that  of 
the  air  front.  Therefore  the  air  probably  does  not  remain  at  the  extremely  high 
temperatures  shown  in  Table  II  for  any  significant  length  of  time,  probably  no 
more  than  the  duration  of  a  few  mean  free  paths.  In  case  1,  had  there  been 
enough  laser  intensity  to  maintain  equilibrium,  the  water  would  have  turned  into 
saturated  steam.  Since  this  result  is  outside  the  two  regimes  of  water  that  were 
analytically  modeled,  it  had  to  be  found  iteratively  by  hand  using  temperature, 
pressure,  density,  and  enthalpy  data  from  the  steam  tables.  The  states  of  the 
water  in  the  three  case  studies,  both  at  the  front  and  in  the  interior,  are  shown  in 
the  Fig.  5  for  comparison. 


TABLE  III 


EQUILIBRIUM  BLAST  MODEL  CASE  STUDIES: 

EQUILIBRIUM  CONDITIONS 

CASE 

1 

2 

3 

Spot  diameter 

m 

0.0032 

0.00066 

0.00017 

Spot  area 

m2 

8.0E-6 

3.4E-7 

2.3E-8 

Incident  intensity 

W/m2 

5.6E+11 

1.3E+13 

2.0E+14 

Energy  density 

MJ/m2 

0.03 

1.0 

14 

INTERIOR  CHARACTERISTICS 

Initial  enthalpy  density  rate 

W/m2 

1.7E+11 

5.3E+12 

1.3E+14 

Final  enthalpy  density  rate 

W/m2 

7.7E+12 

1.5E+13 

2.5E-h14 

Equilibrium  temperature 

oc 

602 

1,747 

17,940 

Laser  intensity  absorbed 

W/m2 

7.5E+12 

9.8E+12 

1.1E+14 

Phase  of  water/air  components 

Sat.SteamVGAS  GAS/GAS 

GAS/GAS 

*  Dryness  fraction 

% 

18.9 

FEASIBILITY 

Input  power  distribution 

Air  front 

% 

13.80 

19.08 

31.86 

Waterfront 

% 

3.05 

6.39 

11.10 

Interior  equilibrium 

% 

1341.97 

74.53 

57.04 

TOTAL  % 

(1358.82) 

100.00 

100.00 

Is  the  laser  power  adequate? 

NO 

YES 

YES 

CONDUCTION  LIMITED  CONDITIONS 

Critical  thickness  (1 -dimensional  stage)  m 

3.6E-7 

2.6E-5 

1.5E-4 

(1 -dimensional  stage) 

s 

5.5E-11 

1.4E-9 

2.9E-9 

CHARACTERISTICS  OF  THE  2-DIMENSIONAL  STAGE 

''2D 

s 

1.3E-7 

8.4E-9 

7.4E-10 

^u 

s 

1.3E-7 

8.4E-9 

2.3E-9 

LASER  PULSE  UTILIZATION 

% 

100 

22.39 

6.02 

It  appears  that  the  laser  energy  must  enter  the  blast  volume  through  the  air 
front.  From  the  air  blast  front,  it  is  transported  through  the  interior  towards  the 
underwater  front.  There  are  a  number  of  possible  transport  mechanisms.  The 
two  main  candidates  are  conduction  and  reradiation.  Conduction  in  this  case  is 
taken  to  include  convection.  Under  the  reradiation  category,  there  are  several 
candidates,  including  simple  scattering,  Raman  scattering,  and  stimulated 
emissions.  We  will  only  consider  the  conduction  mechanism  in  this  study. 

Let  us  now  consider  the  conduction  mechanism  in  greater  detail. 
Conduction  will  take  place  naturally  since  the  temperature  at  the  underwater 
front,  as  the  numerical  results  show,  is  always  significantly  less  than  that  at  the 
air  front.  Furthermore,  any  drop  in  pressure  at  the  underwater  front  will  be 
accompanied  by  a  further  temperature  reduction,  which  will  increase  the 
temperature  gradient.  Therefore,  as  long  as  the  distance  between  the  two  fronts 
is  small  such  that  the  thermal  resistance  is  negligible,  the  process  will  be  a 
dynamic  and  stable  one  in  which  the  input  laser  energy  tends  to  be  distributed 
in  a  way  that  will  maintain  the  equilibrium.  The  thermal  conductivity  k  of  steam 
and  air  may  be  approximated  by  a  linear  temperature  dependence, 

k  =  kj  +  kg  ,  (42) 

where  k,  and  kg  are  constants. 

For  air,  it  was  found  that  using  the  values  (0.000062,  0.0054)  for  (ko,ki)  in 
SI  units,  the  above  conductivity  function  gives  values  that  are  very  close  to  the 
tabulated  values^^  over  a  wide  range  of  conditions.  The  tabulated  values  of  the 
thermal  conductivity  of  water  do  not  increase  monotonically  with  temperature, 
but  within  ±20%  they  can  be  represented  by  a  single  constant  value  of 
0.6  W/m/°K.  The  conductivity  of  steam  at  low  temperatures  has  a  very 
complicated  dependence  on  pressure  and  temperature,  but  curves  constructed 
from  the  tabulated  values  suggest  that  at  higher  temperatures  it  also  converges 
to  a  linear  temperature  dependence  of  the  form  shown  in  Eq.  (42);  it  was  found 
that  values  for  (ko,ki)  of  (0.000062,  0.012)  could  be  expected  to  give  a  good  fit. 


The  comparisons  between  tabulated  values  and  the  analytical  approximations 
are  shown  in  Fig.  6. 


The  necessary  heat  flow  rate  per  unit  area  is  given  by  the  power  density 
required  to  drive  the  underwater  front  and  to  maintain  equilibrium  in  the  blast 
fluids  (Iqu  +  log)-  This  is  related  to  the  temperature  gradient  dT/dx  by 

dT/dx  =  (lou  +  log)/k  .  (43) 

where  x  is  vertical  distance.  Substituting  for  k  from  the  above  equation  and 
integrating,  a  relationship  is  found  between  the  temperatures  Tgo  and  Tyo  at  the 
air  and  underwater  fronts,  the  power  density  (Iqu  +  log),  and  a  critical  thickness 
Xc,  given  by 

(lou  +  log)  Xe  =  l<o(Tao  '  T„„)  +  l<l(Tao^  -  •  (44) 

For  a  layer  of  a  thickness  less  than  Xq,  the  energy  distribution  within  the  layer 
will  remain  in  dynamic  equilibrium.  When  the  layer  thickness  exceeds  Xc,  the 
free  flow  of  energy  from  the  air  front  to  the  underwater  front  will  be  impeded  and 
the  pressure  at  the  underwater  front  will  decline.  With  the  air  and  underwater 
fronts  moving  apart  at  a  velocity  of  (Dgo+Dyo),  the  critical  thickness  Xg  will  be 
reached  at  an  estimated  time  of  after  the  start  of  the  process,  where  is 
simply  given  by 

\  =  V(Dao  +  Duo)  .  (45), 

Computed  values  of  for  the  one-dimensional  stage  are  shown  in  Table  III  for 
the  three  case  studies. 

After  the  critical  thickness  is  reached,  the  energy  distribution  within  the 
blast  volume  will  be  distorted.  An  increasingly  smaller  proportion  of  the  energy 
will  reach  the  underwater  front.  The  pressure  at  the  underwater  front  will  drop 
giving  rise  to  a  pressure  gradient  within  the  blast  volume  which  will  tend  to  drive 
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the  blast  gases  into  the  water.  This  is  a  likely  mechanism  for  the  subsequent 
formation  of  an  underwater  gas  bubble.  The  air  front  may  be  expected  to 
continue  to  travel  at  more  or  less  its  original  speed.  The  underwater  front  will 
slow  down.  Since,  as  the  numerical  results  in  Table  II  show,  the  underwater 
front  velocity  and  temperature  are  small  compared  to  those  of  the  air  front,  it  can 
be  deduced  that  the  layer  thickness  will  be  mainly  determined  by  the  air  front 
velocity  Dgo-  Consequently,  the  power  density  ly  at  the  underwater  front  will 
decline  as  the  reciprocal  of  time.  Thus,  the  power  density  1^  may  be  expressed 
in  terms  of  the  initial  power  density  Iqu  divided  by  a  correction  factor 

>u  =  lou/Pk  :  for  t  >  ,  (46) 

where 

Pk  =  t/tk  .  (47) 

4.  The  Two-Dimensional  Stage 

As  the  two  blast  fronts  move  away  from  each  other,  the  thickness  of 
the  interior  volume  will  increase.  In  addition  to  its  vertical  expansion,  the  blast 
volume  also  expands  laterally  outward.  It  is  reasonable  to  ignore  the  lateral 
expansion  and  treat  the  process  as  one-dimensional  until  the  front 
displacements  become  comparable  to  their  diameters.  As  the  vertical  and 
lateral  displacements  increase,  the  frontal  area  will  also  increase.  This  will 
cause  the  rate  of  mass  absorption  to  be  significantly  increased  and  thus  the 
input  laser  power  will  have  to  be  diluted  over  an  increasingly  larger  mass  input 
rate.  The  idealized  geometry  is  shown  in  Fig.  7,  which  shows  the  air  and 
underwater  fronts  expanding  both  vertically  and  laterally.  The  result  is  a 
reduction  in  the  intensity  available  to  drive  the  undenA^ater  front.  The  effects  on 
the  pressure  at  the  underwater  front  are  of  particular  interest  since  it  will  have  a 
direct  impact  on  the  acoustic  output  of  the  process. 
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To  account  for  the  reduction  in  pressure  at  the  underwater  front,  the  energy 
equation  of  Eq.  (13)  as  it  applies  to  the  underwater  front  is  now  modified  by 
replacing  the  power  input  per  unit  area  term  l^,  with  a  time  dependent  term  l^  for 
the  underwater  blast.  Let  lu  be  expressed  in  terms  of  its  initial  value  Iq^  and  a 
correction  term  by,  giving 

lu/(roD)  =  louAroDbu)  =  w  +  O.Su^  -  O.SD^  .  (48) 

If  it  can  be  assumed  that  the  thermal  resistivity  is  negligible,  then  the  input 
energy  will  be  distributed  in  a  way  that  will  maintain  pressure  and  thermal 
equilibrium  within  the  volume.  Then,  the  parameter  bu  is  given  by  the  ratio 
between  the  mass  absorption  rates  between  the  one-dimensional  stage  at  time 
t=0+  and  at  a  general  time  t.  With  reference  to  the  geometry  shown  in  Fig.  7,  it 
can  be  seen  to  be  given  by 

bu(  DuoW,o+  )7tdo^/4 

a  Wu  (  Durt(do  +  2hu)^/4  +  n  Du(do  +  2hu)hu ) 

+  Wa  (  Da7l(do  +  2ha)^/4  +  7t  Da(do  +  2ha)ha  )  .  (49) 

The  terms  Wg  and  are  weights  related  to  the  relative  energy  consumption 
rates  of  the  fronts  and  the  interior  given  by 

Wa  =  ((IAMIg/'u)(Da/(Da  +  D,)) 

and 

Wu  =  (1  +  (IA)(Du/(Da  +  D J)  .  (50) 

The  terms  Wg,  Wu,  la,  lu.  Ig.  Da-  ^*^^1  Du  are  functions  of  time  t  while  Wgo,  Wuo,  Igo. 
■uo.  lgo>  ^ao.  ^^icl  Duo  3'’®  ^^eir  initial  values  in  the  one-dimensional  stage  at  t=0.^; 
do  is  the  initial  diameter  of  the  frontal  area,  which  is  also  equal  to  the  diameter  of 


the  laser  beam,  and  hg  and  are  the  vertical  distances  traveled  by  the  fronts, 
given  by 

ha  =  Da  dx 
and 

hu  =  f  Du  dx  .  (51) 

0 

Let  us  model  Dg  and  Dy  as  functions  that  are  approximately  constant  in  the  one¬ 
dimensional  stage,  defined  as  the  interval  0<  t  ^X2d  where  X2d  is  taken  to  be  a 
time  constant  characteristic  of  the  duration  of  the  one-dimensional  stage.  For 
t>X2D  let  Dgand  Dy  asymptotically  decay  as  (t/x2D)®.  where  s  is  a  negative  power. 
Using  this  model,  the  time  integral  of  Dg  may  be  approximated  by 

J  Dg  dx  =  tDgo(x2D/t+  (T/X2Df/i  dx) 

=  tDgo(X2D/t  -  (X2d/I)/(1+S)  +  (t/X2D)Ml+S)  )  .  (52) 

and  similarly  for  that  of  Dy.  This  model  will  also  conveniently  make  the  weights 
Wg  and  Wu  time  independent.  Assuming  that  0>s>-1  and  large  t,  the  time 
integral  of  Dg  may  be  approximated  by  the  asymptote 

f  Dg  dx  tDgo(t/X2D)"/(1+s)  .  (53) 

0 

A  similar  result  is  also  obtained  for  the  time  average  of  Du-  Substituting  the 
asymptotic  forms  of  Dg  and  Dy  and  their  integrals  into  Eq.  (49),  the  asymptotic 
form  of  Pu  is  given  by 


p„  ^  W„„D„„3)(t/t2Dr*^^2on1+S) 


X  (  DuoW^  +  W^Dao  )  'd„'^  .  (54) 

Using  Iqu/Pu  instead  of  Iq  in  the  expression  for  front  velocity  of  Eq.  (1),  it  is 
deduced  that  varies  as  Since  the  above  equation  shows  that  pu  is 

asymptotically  proportional  to  t  ,  equating  the  powers  of  t,  it  is  deduced  that 

s  =  (-1/3)(3s+2)  .  (55) 

Solving,  it  is  found  that  the  value  of  s  is  -1/3,  which  justifies  our  initial 
assumption  of  its  range  of  values.  Substituting  for  s  in  Eq.  (54),  it  is  found  that 
Pu  must  vary  linearly  with  t.  Then,  using  Iqu^Pu  instead  of  Iq  in  Eq.  (15),  it  is 
deduced  that  the  pressure  Pu  at  the  underwater  front  must  vary  as 

Pu  ^  Puo('C2D/t)^  :  fort>t2o  ,  (56) 

where  p^g  is  the  initial  pressure  in  the  one-dimensional  stage. 

To  obtain  an  approximate  estimate  of  T2d  *^5  equate  the  asymptotic  form 
of  Pu  to  the  initial  value  of  Pu  which  is  one,  replace  t  with  T2d.  and  solve,  giving 

X2D  =  do(1  +S)[(DuoWuo  +  WaoDao)/(12(WaoD3u3+  Wu^Duo^))]’''^  ■  (57) 

Calculated  values  of  X2D  three  case  studies  are  shown  in  Table  III. 


Let  us  now  consider  the  case  where  the  blast  enters  the  two- 
dimensional  stage  and  then  enters  the  conduction  limited  condition,  i.e.,  the 
case  where  X2D<T^^  for  time  t>Xk.  This  is  applicable  when  the  laser  beam  is 
narrow;  thus  the  blast  is  apt  to  enter  the  two-dimensional  stage  before  its 
dimensions  approach  the  critical  thickness  of  the  conduction  limited  condition. 


i 


The  correction  term  p^'  for  the  conduction  limited  condition  will  differ  by  a 
constant  factor  from  p^  as  given  by  Eq.  (47)  because  the  decline  in  pressure  due 
to  the  two-dimensional  expansion  will  alter  the  time  constant  of  the  conduction 
limited  condition.  Otherwise,  the  intensity  reductions  due  to  expansion  and 
thermal  conduction  limitation  are  expected  to  be  independent,  and  hence 
multiplicative.  Using  the  asymptotic  forms  of  p^  and  p^',  the  power  density  at  the 
underwater  front  1^  will  be  given  by 

lu  =  lou/(PuPk’) 

louM)^  :  fort»Xu  ,  (58) 

where  is  the  new  time  constant  for  the  combined  conduction  and  expansion 
effects.  Replacing  Iq  in  Eq.  (15)  with  l^  from  the  above  equation,  it  is  seen  that 
the  pressure  will  decline  as  Thus,  the  underwater  front  pressure  p^  is 
asymptotically  given  as 

Pu  ^  Puo  :  t « tu 


PuoM) 


:  for  t»Xu 


To  obtain  an  estimate  of  x^  it  is  necessary  to  return  to  Eq.  (44)  and  replace  the 
constant  terms  (lou+  log).  Tgo,  and  T^o  with  their  time  dependent  counterparts  (1^+ 
Ig),  Tg,  and  T^.  From  the  above  analysis  of  (Iqu-i-  l^g),  Tgo,  and  T^q,  it  is  deduced 
that  their  asymptotic  forms  are 


(lu+  Ig)  (lou+  log)('^2D'^) 


Ta  Tgo  (X2D/t)^ 

1’uo  ('^20^)^'^ 


;  fort>X2D 

:  for  t  >  X2d 
:  fort>X2D 


Substituting  x^,  for  t  and  replacing  ((0^,+  l^g),  Tg^,  and  T^^^  with  (1^+  Ig),  Tg. 
and  Tu,  in  Eq.  (44),  a  relationship  between  x,.  and  is  found. 


.v?.  '.-■vivi  .-ivv’.v;' .-VI' 
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(Iqu  log)(^2D^"^u) 


=  KiC^ao  * '^uo)(^2D^^u)  '*’  ^iC^ao  ‘  ^uo  )(^2D^^u) 


A  second  relationship  between  x^,  and  is  given  by 


f'^u 

Xc  =  J  Da  +  Dudx 

0 

2/3 

(Dao  +  Duo)t2D(1  +  3((xyt2D)  '  1)/2)  i  for  T,  >  X2D  •  (64 

Eliminating  x^  from  these  two  equations,  a  cubic  equation  in  {x2D/tu)^^  is  found, 


(Iqu  iog)'^2D{^ao  ^uo)(  3  *  ('^2D^^u) 

=  ko(Tao  * 'ruo)('C2D'^^u)  ^iC^ao  *  ^uo  )i'^2D^%)  ^3,  (65) 

which  may  be  solved  using  standard  solutions. 

In  the  case  where  the  blast  enters  the  conduction  limited  condition  while 
still  in  its  one-dimensional  stage,  i.e.,  X2Di'rk>  underwater  front  pressure  will 
decline  according  as  until  the  two-dimensional  stage  is  reached.  It  will 
reach  the  two-dimensional  stage  at  a  time  not  significantly  different  from  X2d  as 
defined  in  Eq.  (57).  This  is  because  the  arrival  time  of  the  two-dimensional 
stage  is  mainly  controlled  by  the  velocity  of  the  air  blast  front  since  it  is  very 
much  faster  than  the  underwater  front;  and,  since  the  air  front  velocity  is  largely 
unaffected  by  the  conduction  limited  condition,  it  is  expected  that  X2d  'S  also 
largely  unaffected.  Therefore,  in  this  case,  the  time  constant  for  the  combined 
conduction  and  expansion  effects  is  simply  equal  to  X2d- 


The  values  of  for  the  three  case  studies  are  shown  in  Table  III.  It  is  seen 
that  the  values  of  x^,  can  be  very  small,  ranging  from  0.13  ^s  in  case  1  to  2.3  ns 
in  case  3. 

The  above  results  suggest  that  there  is  a  definite  time  window  in  which  the 
laser  energy  will  have  the  maximum  effect.  This  time  window  is  expected  to  be 
closely  related  to  By  integrating  the  pressure  Pu  as  defined  in  Eq.  (59),  over 
all  time,  a  simple  expression  for  the  total  impulse  delivered  to  the  water  by  the 
blast  is  obtained. 

Impulse  =  J  Pg  dt 
0 

=  4tgPuo  .  (66) 

It  can  be  shown  that  more  than  half  of  the  impulse  is  delivered  between  t=0  and 
t=4Tu.  Therefore,  the  usefulness  of  any  laser  energy  that  arrives  after  4Tg  will  be 
much  reduced.  Thus,  dividing  the  laser  pulse  length  by  4Tg,  an  estimate  of  the 
percentage  of  the  laser  energy  that  is  most  useful  may  be  estimated.  The 
results  for  the  case  studies  as  shown  in  Table  III  indicate  that,  while  in  case  1  all 
the  laser  energy  is  useful,  in  case  3  the  useful  proportion  is  only  6.02%. 

6.  Blast  Evolution 

From  the  above  analysis  and  numerical  case  studies,  it  is  possible  to 
piece  together  a  qualitative  picture  of  the  evolution  of  the  blast  process. 

In  the  one-dimensional  stage,  two  fronts  will  form  at  the  interface  and  begin 
to  move  apart.  With  its  much  higher  equilibrium  velocity,  the  air  blast  front  will 
move  away  from  the  interface  at  a  faster  rate  than  the  underwater  front.  Thus, 
the  blast  volume  will  be  mainly  above  the  water  surface,  as  illustrated  in 
Fig.  8(a).  The  blast  will  send  a  shock  front  into  the  water.  At  X2d  the  blast 
enters  the  two-dimensional  stage,  in  which  dilution  of  the  input  power  from  the 
laser  over  an  increasingly  larger  frontal  area  and  interior  volume  becomes 
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significant.  Therefore,  the  fronts  will  continue  to  expand  but  at  increasingly 
slower  rates  as  the  pressure  declines.  The  volume  evolves  from  an  initial 
pancake  shape  into  a  more  rounded  shape,  as  illustrated  in  Fig.  8(b).  The 
increase  in  the  thermal  resistance  due  to  the  widening  gap  will  become 
significant  after  time  and  result  in  an  increasingly  smaller  proportion  of  the 
laser  energy  delivered  to  the  undenwater  front.  The  part  of  the  air  front  in  direct 
contact  with  the  laser  beam  will  continue  to  travel  at  a  high  speed  while  the 
remaining  fronts  lose  pressure  and  slow  down,  as  illustrated  in  Fig.  8(c).  The 
pressure  gradient  which  subsequently  develops  will  tend  to  drive  gas  into  the 
water.  The  gas  driven  into  the  water  by  the  pressure  gradient  will  collect  into  a 
bubble  which  may  oscillate,  causing  one  or  more  subsequent  bubble  pulses. 

High  speed  photography  of  the  interaction  of  a  high  power  laser  with  a 
water  surface  by  Emmony,  Geerken,  and  Straaijer®  cleariy  shows  the 
development  of  the  later  stages  of  the  blast.  Their  photographs  show  the 
development  of  a  spherical  blast  whose  center  is  well  above  the  interface. 
Later,  when  the  effects  of  thermal  resistivity  become  significant,  they  show  the 
part  of  the  air  front  that  is  in  direct  contact  with  the  laser  beam  continuing  to  rise 
to  great  heights,  while  the  expansion  rate  of  the  main  volume  declines.  Finaliy, 
a  gas  bubble  is  driven  below  the  water  surface.  Their  photographs  also  show 
subsequent  bubble  activity  which  is  beyond  the  scope  of  this  blast  model. 

7.  Underwater  Acoustic  Cutout 

In  the  case  where  the  blast  is  weak,  the  velocity  of  the  underwater 
front  will  be  approximately  equal  to  that  of  small  amplitude  acoustic  waves.  The 
effective  intensity  1^  will  drop  below  the  minimum  requirement  early  in  the  one¬ 
dimensional  stage.  At  this  point,  the  front  will  have  moved  a  negligible  distance, 
and  therefore  all  the  significant  pressure  changes  may  be  considered  as  taking 
place  at  the  interface.  As  an  acoustic  source,  it  may  be  modeled  as  a  planar 
pressure  transducer.  Since  the  pressure  is  weak,  linear  superposition  may  be 
applied,  and  thus  the  acoustic  signal  at  a  distant  point  may  be  simply 


constructed  through  a  Green's  function  integral.  Thus  the  acoustic  pressure  p(r) 
at  a  point  r  is  given  by  the  surface  integral 


P(')  =  _/”_/“(<iPu(t-|r-s|/Coy<«V(|r-s|Po)  ds  ,  (67) 

where  Pu(t  -  |r‘S|/Co )  is  the  time  retarded  underwater  pressure  and  s  is  a  general 
point  on  the  interface.  This  approach  was  taken  in  the  previous  study.22  This 
allowed  the  optoacoustic  process  to  be  studied  without  the  complication  of 
nonlinear  acoustic  propagation  effects. 

In  the  case  of  a  strong  blast,  the  underwater  front  will  continue  to  be  driven 
well  after  1^,  that  is,  well  into  the  conduction  limited  condition.  To  obtain  a  better 
understanding  of  the  acoustic  generation  process,  let  us  look  more  closely  at 
the  development  of  the  undenivater  front,  particularly  at  the  heat  flow  to  the  front 
and  the  shock  front  itself. 

Let  us  consider  the  blast  gases  behind  the  underwater  front.  In  the 
conduction  limited  condition,  the  energy  deposited  by  the  laser  beam  cannot  be 
instantaneously  distributed.  Therefore  a  significant  concentration  of  heat 
energy  will  build  up  in  the  collection  volume.  Since  energy  collection  takes 
place  at  the  air  front,  the  collection  volume  is  the  volume  occupied  by  the  laser 
beam  that  has  been  swept  by  the  air  front.  For  a  circular  laser  beam  directed 
vertically  downward,  the  collection  volume  is  a  vertical  cylinder  starting  at  the 
interface  and  terminating  at  the  air  front,  as  shown  in  Fig.  9.  Therefore  the 
governing  process  controlling  the  development  of  the  underwater  front  Is  the 
heat  diffusion  from  the  cylindrical  reservoir,  particularly  through  its  bottom  end. 
The  problem  of  heat  diffusion  is  one  that  satisfies  Laplace's  equation,  and  may 
be  solved  by  standard  methods  in  the  static  case.  From  the  geometry  of  the 
problem,  it  can  be  seen  that,  for  a  point  on  the  underwater  front  directly  below 
the  reservoir,  the  downward  diffusion  from  the  curved  sides  of  the  reservoir  will 
be  small  compared  to  that  from  the  flat  bottom  end.  Neglecting  the  diffusion 
from  the  curved  sides,  the  problem  becomes  one  of  diffusion  from  a  circular 
plate.  In  this  case,  the  solution  of  the  static  problem  is  very  simple:  Using  a 
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conformal  transformation  technique,  it  may  be  shown  that  the  lines  of  flux 
are  confocal  hyperbolae  that  are  orthogonal  to  confocal  ellipsoid  constant 
temperature  surfaces,  as  illustrated  in  Fig.  9.  Thus  the  underwater  constant 
temperature  surfaces  will  tend  to  conform  to  the  confocal  ellipsoid  surfaces, 
since  any  deviation  from  this  basic  shape  would  be  accompanied  by  a 
corrective  change  in  the  heat  input  rate.  The  foci  are  located  on  the  rim  of  the 
circular  heat  reservoir.  Assuming  that  the  diameters  of  the  heat  reservoir  and 
the  laser  beam  are  approximately  equal,  the  general  equation  of  the  constant 
temperature  surface,  for  a  laser  beam  of  diameter  do,  is  given  by 

z2  =  (q/2)2  -  (do/2)2  -  r2(  1  -  (do/q)2)  ,  (68) 

where  z  is  the  depth  of  a  point  on  the  constar^t  temperature  surface  and  r  its 
radial  distance  from  the  vertical  axis  of  symmetry:  q  is  the  independent 
parameter. 

On  the  vertical  axis  of  symmetry,  i.e.,  at  r=0,  the  radius  of  curvature  Rp  of 
the  ellipsoid  is  given  by 

=  -1/(d2z/dr2) 

=  (q/2)  ( 1  -  (dVq)2)-''2  .  (69) 

Differentiating  Rg  with  respect  to  q  and  setting  the  result  equal  to  zero,  it  is  found 
that  the  value  of  Rc  has  a  minimum  of  dg,  which  occurs  at  q  =  doV2.  Substituting 
for  q  in  Eq.  (68)  and  setting  r  to  zero,  it  is  found  that  the  minimum  radius  of 
curvature  occurs  at  z=(iJ2. 

Initially,  the  underwater  front  is  modeled  as  a  perfect  surface  of 
discontinuity,  where  the  fluids  enter  and  exit  an  infinitely  thin  transformation 
zone.  In  the  one-dimensional  stage,  the  underwater  front  is  a  shock  front, 
driven  through  the  water  by  hot  gases  from  behind,  traveling  at  a  constant 
speed  and  in  a  state  of  dynamic  equilibrium.  As  the  front  enters  the  conduction 


limited  condition,  its  input  intensity  1^  declines  and  therefore  it  is  no  longer  able 
to  process  the  incoming  fluids  at  the  same  rate.  As  a  consequence  two  things 
will  happen;  (1)  the  front  will  slow  down  and  (2)  it  will  "thicken",  i.e.,  the  exit 
surface,  or  gas  boundary,  will  begin  to  lag  behind  the  shock  front,  eventually 
giving  way  to  the  complete  separation  of  a  self-contained  shock  pulse  from  the 
gas  bubble.  There  will  be  strong  pressure  and  particle  velocity  gradients  in  this 
zone.  Since  the  particle  velocities  in  a  strong  shock  front  are  quite  substantial, 
there  is  a  case  for  treating  the  dynamics  in  this  zone  in  terms  of  a  flow  field 
rather  than  a  wave  field.  Numerical  results  of  tj  from  the  case  studies  shown  in 
Table  II  indicate  that  the  change  in  the  density  of  water  caused  by  the 
underwater  shock  front  is  extremely  small.  Therefore  the  velocity  field  of  the 
particle  motion  may  be  regarded  as  nearly  incompressible.  Since  the  flow  field 
of  an  inviscid  and  incompressible  fluid  also  obeys  Laplace's  equation,  it  is 
expected  that  the  evolution  of  the  shock  front  will  also  have  a  tendency  to 
conform  to  the  confocal  ellipsoid  model,  as  illustrated  in  Fig.  9.  In  the  limit  as 
the  parameter  q  in  Eq.  (68)  tends  to  infinity,  the  ellipsoid  will  become  a  sphere. 
Thus,  the  ellipsoid  model  will  allow  the  shock  front  to  smoothly  develop  into  a 
spherical  front  in  the  farfield. 

Therefore,  it  is  argued  that  the  ellipsoid  model  is  the  asymptotic  model  that 
the  underwater  front  will  tend  towards,  in  the  limiting  case  of  extremely  high 
heat  flow  rates  in  the  blast  volume  and  high  particle  velocities  behind  the  shock 
front.  In  practice,  it  serves  as  a  useful  but  approximate  model. 

From  studies  of  underwater  explosions,^^  there  are  only  a  few  main 
characteristic  quantities  that  define  the  final  spherically  spreading  shock  front, 
the  initial  peak  pressure,  the  radius  of  curvature,  and  the  rate  of  decay  (or  time 
constant).  The  theoretical  results  developed  in  the  foregoing  section  may  be 
used  to  estimate  these  characteristic  quantities.  The  propagation  of  the  shock 
wave  may  then  be  modeled  as  that  of  a  simple  spherical  shock  for  which  the 
pressure  and  time  constant  variations  have  been  theoretically  analysed  by 
several  researchers.  For  the  purposes  of  this  analysis,  the  results  of  Rogers^'' 
will  be  used.  The  expressions  for  shock  peak  pressure  p^  and  time  constant  x, 


in  terms  of  the  initial  values  of  range  Rq.  pressure  Po,  and  time  constant  Tq,  are 
given  by 


Pm  =  Po{[1 +2(Ro/Xo)ln(R/Ro)]''^-1}/[(R/Xo)ln(R/Ro)]  .  (70) 

t  =  T,  [1  +  2(R^XR/Ro)]''^  .  (71) 

where  the  characteristic  length  Xq  is  given  in  terms  of  the  coefficient  of 
nonlinearity  p  (value  3.5  for  water)  and  the  initial  values  by 

=  PoCo^VPPo  •  (72) 

The  initial  peak  pressure  Po  will  be  taken  as  equal  to  the  initial  blast 
pressure  as  formulated  in  Eq.  (29).  The  initial  range  Rq  will  be  taken  as  equal  to 
the  radius  of  curvature  Rc  of  the  ellipsoid  model  at  its  minimum  point.  In  the 
vertically  downwards  direction,  it  was  found  through  Eq.  (69)  to  be  equal  to  the 
diameter  of  the  laser  beam  d^,  giving 

Ro  =  do  ■  (73) 

The  initial  time  constant  will  be  obtained  from  the  asymptotic  decay  rates  as 
given  by  Eqs.  (59)  and  (66)  by  equating  the  impulses  of  the  exponential  pulse  of 
the  Kirkwood-Bethe  model^^  and  that  of  the  above  model.  The  impulse  Ikb  in 
the  Kirkwood-Bethe  model  is  given  by 

foo  Too  W- 

Ikb  =  J  Pu  dt  =  J  p^o  e  «  dt  =  p^oXo  •  (74) 

The  impulse  Ig  in  the  above  model  is  given  by 

-4/3 

Ig  =  J  Pu  dt  PuoXu  +  -  J  Puo(t/^u)  dt  =  4puotu  .  (75) 

0 


Equating  Ikb  and  Ig,  the  result  is  simply 


Applying  the  commonly  used  approximation  that  the  final  shock  pulse  is 
also  similar  to  an  exponential  function,  the  values  of  the  peak  pressure  p^  and 
time  constant  x  may  be  used  to  obtain  approximate  but  simple  estimates  of  the 
magnitude  and  the  roll-off  of  the  shock  pulse  energy  density  spectrum  |P(a))|^, 
giving 

|P(w)l^  ==  (Pm'c)^/(  1  +  ) :  forO«a)  .  (77) 

Numerical  results  for  the  three  case  studies  of  the  initial  and  final  values  of 
the  key  parameters  are  shown  in  Table  IV;  in  the  example  shown,  the  final 
values  refer  to  a  point  of  observation  1  m  directly  below  the  blast.  The  figures 
for  case  1  are  in  brackets  since  they  would  only  apply  if  the  laser  intensity  in 
case  1  were  adequate,  i.e.,  about  14  times  (1359%)  higher  than  the  stated 
value  in  Table  III.  It  is  interesting  to  note  that,  while  the  initial  peak  pressure  is 
highest  in  case  3,  the  final  pressure  is  not.  This  is  due  to  its  smaller  initial  radius 
of  curvature.  Of  the  three  cases,  the  case  that  is  predicted  to  deliver  the  highest 
peak  pressure  and  spectrum  density  level  is  case  2.  Case  1  proved  to  be  a 
non-starter;  the  case  1  values  in  brackets  are  not  applicable  because  they  are 
computed  for  a  laser  power  that  is  1 4  times  higher  than  the  stated  value. 

In  Fig.  10  a  plot  of  the  predicted  acoustic  peak  pressure,  for  the  point  1  m 
directly  below  the  blast,  as  a  function  of  intensity  is  shown.  The  laser  power  is 
assumed  to  be  constant  and  the  intensity  is  varied  by  changing  the  spot  size. 
Two  curves  are  shown,  one  for  a  laser  power  of  4.5  MW  which  is  the  power 
delivered  by  our  Nd:glass  laser,  and  a  second  curve  for  a  laser  power  that  is  ten 
times  higher.  It  was  found  that  above  an  intensity  of  300  TW/m^  the  fluid  exiting 
the  underwater  front  is  predicted  to  be  steam.  Below  this  intensity,  it  is  liquid 
water  which  is  subsequently  turned  into  steam  on  contact  with  the  hot  air  from 
the  air  front.  It  is  seen  that  below  intensities  of  about  10  TW/m^,  there  is 
insufficient  intensity  to  maintain  the  equilibrium  within  the  blast,  that  is,  to  drive 


TABLE  IV 


EQUILIBRIUM  BLAST  MODEL  CASE  STUDIES:  ACOUSTIC  OUTPUT 


CASE  STUDIES 

1 

2 

3 

Spot  diameter 

m 

0.0032 

0.00066 

0.00017 

Spot  area 

m2 

8.0E-6 

3.4E-7 

2.3E-8 

Incident  intensity 

W/m2 

5.6E+11 

1.3E+13 

2.0E+14 

Energy  density 

MJ/m2 

0.03 

1.0 

14 

ACOUSTIC  OUTPUT : 

^0 

s 

{5.4E-07) 

3.4E-08 

9.0E-09 

Initial  peak  pressure 

Pa 

(1.3E+7) 

1.3E-t-8 

1.1  E+9 

*0 

m 

(4.3E-02) 

2.7E-04 

8.3E-06 

Rq  (vertical  case) 

m 

(3.2E-3) 

6.6E-4 

1.7E-4 

R 

m 

(1) 

1 

1 

Time  constant  at  R;  t 

s 

(7.3E-7) 

2.1E-7 

1.7E-7 

Peak  pressure  at  R:  Pm 

Pa 

(3.5E+4) 

2.4E+4 

1.9E+4 

Spectrum  level  dB  re  1  W/m2/Hz  re  1  m 

(-94) 

-108 

-112 

Upper  cutoff  frequency 

MHz 

(0.21) 

0.78 

0.94 

The  figures  for  case  1  are  in  brackets  because  they  represent  the  values  that  would  have  been 
obtained  had  there  been  enough  laser  power  to  drive  the  process. 
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the  fronts  at  the  maximum  entropy  condition  and  sustain  the  pressure 
equilibrium  within  the  blast  interior.  It  is  possible  that  a  pressure  wave  may  still 
be  generated  outside  blast  equilibrium  conditions;  therefore,  in  practice,  the 
pressure  is  not  expected  to  decay  quite  so  abruptly  with  intensity.  It  is  seen  in 
Fig.  10  that  acoustic  pressure  maxima  occur  at  the  minimum  intensity  (maximum 
spot  size)  required  to  sustain  an  equilibrium  blast. 

In  Fig.1 1 ,  a  similar  plot  of  the  low  frequency  spectrum  levels  is  shown.  The 
low  frequency  spectrum  estimate  is  based  on  the  assumption  of  an  exponential 
shock  function.  Under  this  assumption,  the  low  frequency  end  of  the  acoustic 
pressure  Fourier  transform  is  approximately  constant,  given  by  the  product  of 
the  peak  pressure  and  the  time  constant. 

Assuming  that  the  radiated  acoustic  energy  approximately  follows  a  dipole 
directivity  pattern,  the  total  acoustic  energy  radiated,  at  a  range  of  1  m,  was 
estimated.  The  acoustic  energy  was  divided  by  the  useful  laser  energy,  giving  a 
measure  of  the  energy  conversion  efficiency  for  a  range  of  1  m;  the  useful  laser 
energy  is  defined  as  the  stated  laser  power  times  the  useful  pulse  length  which 
is  estimated  to  be  81^,  where  is  the  pressure  decay  time  constant  as  defined 
by  Eqs.  (65)  and  (66).  The  plots  of  efficiency  as  a  function  of  spot  size  and 
intensity  are  shown  in  Fig.  12(a)  and  (b).  Also  shown  are  plots  of  efficiency  for 
the  acoustic  energy  in  the  band  0-50  kHz;  this  band  was  chosen  because  it 
covers  a  broad  range  of  sonar  applications.  It  is  seen  in  Fig.  12(a)  that,  below 
the  minimum  spot  size  required  to  generate  steam  within  the  blast  interior,  the 
fullband  efficiencies  are  practically  independent  of  the  laser  power,  with  a 
minimum  value  of  approximately  0.04%.  It  is  also  seen  that  the  0-50  kHz  band 
efficiencies  are  about  10  times  lower  and  increase  only  slowly  with  laser  power. 
It  is  also  seen  in  Fig.  12(b)  that  the  efficiency  in  the  0-50  kHz  band  is 
independent  of  laser  power  above  intensities  of  about  100  TW/m^. 
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III.  EXPERIMENT 


A.  Apparatus 

The  apparatus  is  shown  in  Fig.  13.  The  laser  used  is  an  Apollo  Nd:glass 
laser  model  22/A  with  a  maximum  pulse  energy  of  22  J.  The  laser  beam  was 
directed  vertically  down  by  a  prism  and  focused  on  the  water  surface.  With 
appropriate  topside  lenses,  it  is  capable  of  producing  a  spot  size  of 
approximately  8  mm^  (diameter  3.2  mm)  at  the  water  surface.  With  converging 
lenses  nearer  the  water  surface,  spot  sizes  of  0.023  and  0.342  mm^  (diameters 
0.17  and  0.66  mm)  were  also  obtained.  A  short  laser  pulse  was  achieved  by 
using  a  Inrad  model  2-016  Q-switch  driver  to  modulate  the  Pockels  cell  in  the 
laser  cavity.  The  laser  pulse  was  approximately  0.15  ps  long  and  a  train  of  up 
to  five  pulses  at  a  maximum  repetition  rate  of  10,000  pps  could  be  delivered. 
The  maximum  laser  pulse  energy  delivered  at  the  water  surface  per  pulse  was 
0.7  J.  Silicon  PIN  photodiodes  were  used  to  monitor  the  laser  pulse  through 
fiber-optic  probes.  The  acoustic  detectors  used  included  a  E8-59  hydrophone 
and  a  PCB138A01  tourmaline  shock  probe.  The  acoustic  signal  was  recorded 
at  a  number  of  different  observation  angles  q.  The  control  circuit  determined  the 
number  of  pulses  and  the  timing  between  pulses.  The  data,  both  optical  and 
acoustic,  were  recorded  on  a  Nicolet  4094A  digital  oscilloscope. 

Direct  viewing  of  the  reaction  at  the  water  surface  was  not  possible  due  to 
the  speed  of  the  process  and  the  strong  optical  emission.  Therefore,  snapshots 
of  the  activity  on  the  water  surface  were  made  photographically.  The  apparatus 
is  shown  in  Fig.  14.  The  real  image  of  a  grid  pattern  is  projected  on  the  water 
surface  and  captured  on  film  using  a  Polaroid  camera.  Any  disturbance  on  the 
water  surface  is  manifested  as  a  distortion  in  the  grid  image.  By  strobing  the 
grid  with  a  flash  lamp,  the  disturbances  on  the  water  surface  at  a  particular 
instant  may  be  captured.  Since  the  camera  shutter  was  opened  before  the  start 
of  the  optoacoustic  process,  any  optical  emissions  in  the  visible  range  from  the 
surface  produced  at  the  time  of  the  laser  pulse  arrival  were  also  captured  on  the 
film. 
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Let  us  begin  by  comparing  the  experimental  results  with  the  linear 
thermoacoustic  theory,  which  assumes  that  the  acoustic  signal  is  produced  by 
simple  thermal  expansion  of  the  medium,  since  this  is  the  best  understood 
mechanism  of  optoacoustic  sound  generation.  From  linear  thermoacoustic 
theory,^ the  acoustic  signal  spectral  density  |P(cd)|^  is  given  by 

lP(co)|2  =  |W((o)|2  (kCaDa/ro)2exp(-2aro)  ,  (78) 

where 

W{(o)  =  the  Fourier  transform  of  the  laser  power  function, 
k  =  acoustic  wave  number, 

a  =  acoustic  absorption  coefficient, 

ro  =  range  from  optoacoustic  source, 

Ca  =  thermoacoustic  conversion  constant,  determined  by  physical 
constants,  including  specific  heat,  thermal  expansion  coefficient,  etc.,  and 
□a  =  diffraction  loss,  determined  by  the  geometry  and  orientation  of 

laser  beam,  the  laser  beamwidth,  and  absorption  length. 

Dividing  both  sides  of  Eq.  (78)  by  the  square  of  the  laser  energy,  which  is 
identically  W(0),  we  get 

|P(co)/W(0)(2  =  |W(co)/W(0)|2  (kCaDa/ro)2exp(-2arJ  .  (79) 

The  right  hand  side  of  Eq.  (79)  is  now  independent  of  the  laser  energy 
W(0).  It  is  dependent  on  the  width  of  the  laser  beam  through  the  term  Dg  and  on 
the  shape  of  the  laser  pulse  through  the  normalized  spectrum  function 
I W((o)/W(0)|.  In  the  farfield  and  on  the  vertical  axis  directly  below  the 
optoacoustic  source,  it  can  be  shown  that  the  term  Dg  is  independent  of  the 
laser  beamwidth.  The  shape  of  the  laser  pulse  was  kept  constant  throughout 
the  experiment.  Under  these  conditions,  |P(w)/W(0)|^  is  only  dependent  on 


physical  constants;  therefore  it  makes  a  good  metric  for  studying  the 
optoacoustic  conversion  characteristics. 

In  Fig.  15,  representative  experimental  measurements  of  |P(co)/W(0)|^  are 
shown  along  with  the  theoretical  thermoacoustic  characteristic  curve.  The 
acoustic  sensor  was  at  a  depth  of  1  m  directly  below  the  source  point.  Results 
are  shown  from  four  very  different  laser  surface  energy  density  levels,  i.e.,  laser 
energy  delivered  to  the  water  surface  divided  by  beam  cross-sectional  area  at 
the  water  surface.  The  results  show  that  at  0.013  MJ/m^  there  is  good 
agreement  with  the  thermoacoustic  theory  for  the  most  part.  The  discrepancies, 
wherever  they  occur,  tended  to  exceed  the  theoretical  estimate;  therefore  it  is 
deduced  that  the  linear  thermoacoustic  theory  must  be  close  to  the  limits  of  its 
validity.  At  0.03  MJ/m^  and  higher  energy  densities,  it  is  seen  that  the  acoustic 
outputs  were  significantly  higher  than  that  predicted  by  the  thermoacoustic 
theory.  At  14  MJ/m^,  the  acoustic  output  exceeded  that  predicted  by  the 
thermoacoustic  theory  by  more  than  50  dB. 

At  the  lowest  energy  densities,  the  acoustic  pulse  (not  shown)  was  found  to 
have  the  characteristic  N-shape  of  the  bandlimited  time  derivative  of  a  delta 
function,  which  is  as  predicted  by  the  linear  thermoacoustic  theory.  The 
corresponding  acoustic  pulses  for  the  remaining  energy  density  levels  are 
shown  in  Figs.  16,  17,  and  18.  At  0.03  MJ/m^  and  using  the  E8-59  hydrophone, 
the  pulse  is  a  corrupted  figure  of  N,  as  shown  in  Fig.  16,  indicative  of 
nonlinearities.  At  the  higher  laser  energy  densities,  the  signal  was  considerably 
stronger;  thus  it  was  possible  to  use  the  tourmaline  gauge  shock  probe  and  take 
advantage  of  its  greater  bandwidth.  At  1  MJ/m^  and  14  MJ/m^,  as  shown  in 
Figs.  17  and  18,  the  positive  shock  peaks  are  estimated  to  have  time  constants 
of  about  2  |is.  This  is  probably  due  to  the  limited  response  time  of  the  shock 
probe,  which  is  claimed  by  the  manufacturer  to  be  1 .5  ^is.  Therefore,  the  true 
acoustic  peak  pressure  is  likely  to  be  higher  than  that  shown.  With  hindsight,  a 
faster  probe,  such  as  a  PVDF  film,  would  have  been  a  better  choice.  At 
14  MJ/m^,  as  shown  in  Fig.  18,  there  is  significant  clipping  and  stretching  of  the 
rarefaction  phase  due  to  cavitation  limitation. 
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From  photographic  observations  of  the  water  surface,  obtained  using  the 
apparatus  described  in  Fig.  14,  it  was  found  that  at  energy  densities  0.013  and 
0.03  MJ/m^  there  were  no  detectable  surface  deformations.  An  example  of  the 
photograph  of  the  grid  image  on  the  water  surface  under  these  conditions  is 
shown  in  the  Fig.  19.  As  the  energy  density  levels  approached  1  MJ/m^,  some 
disturbances  were  observed.  Occasionally,  as  shown  in  Fig.  20,  the 
photograph  would  show  secondary  emission  from  the  laser  spot;  the  laser  itself 
is  invisible  to  the  film.  Other  times,  as  shown  in  Fig.  21 ,  there  would  be  a  bright 
spot  and  weak  ripples.  At  1  MJ/m^,  for  most  of  the  time,  the  observed  surface 
image  is  as  shown  in  Fig.  22.  This  photograph  shows  a  bright  and  explosive 
optical  emission  from  the  water  surface  where  the  laser  beam  made  contact.  It 
is  also  seen  that  there  are  several  concentric  ripples.  The  photograph  clearly 
shows  that  an  explosive  reaction  has  taken  place.  At  14  MJ/m^  the  ripple 
structure  is  quite  different,  having  only  one  large  ripple  and  a  plume  rising 
vertically  out  of  the  water,  as  shown  in  Fig.  23.  The  plume  is  not  easily  seen,  but 
its  shadow  is  clearly  visible  in  the  foreground.  In  all  the  photographs,  the  grid 
was  strobed  at  10  ms  after  laser  contact;  therefore  the  surface  ripples  appear  as 
they  were  at  that  instant.  The  optical  emissions,  however,  could  have  occurred 
at  any  time  up  to  30  ms  after  the  laser  arrival  since  the  camera  shutter  was  left 
open  during  that  interval. 

The  many  concentric  ripples  seen  in  Fig.  22  are  of  particular  interest.  A 
mathematical  analysis  of  the  water  waves  produced  by  explosions  was  made 
by  Kranzer  and  Keller.^^  Their  analysis  was  done  in  terms  of  the  impulse 
distribution  on  the  water  surface.  For  the  case  where  the  impulse  is  a  nonzero 
constant  within  a  circular  area  and  zero  outside,  they  showed  that  there  will  be 
a  large  number  of  concentric  ripples  of  approximately  equal  amplitude.  It  has 
been  shown  in  Eq.  (66)  of  the  blast  model  that  most  of  the  impulse  is  delivered 
within  a  time  of  4tu.  In  many  cases,  the  lateral  boundaries  of  the  air  and 
underwater  fronts  will  not  have  significantly  expanded  in  that  time;  from  the 
figures  under  case  2  in  Table  II  and  III,  it  can  be  seen  that  this  is  particularly 
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applicable  to  the  data  at  1  MJ/m^.  Thus,  the  impulse  function  may  be 
approximated  by  a  nonzero  constant  within  a  circular  area  and  zero  outside  it. 
The  photograph  of  Fig.  22  shows  several  concentric  ripples  of  roughly  equal 
amplitude  which  seems  to  be  in  agreement  with  the  mathematical  results  of 
Kranzer  and  Keller. 

Apparently,  this  is  not  the  case  for  the  data  at  14  MJ/m^.  In  this  case,  the 
blast  diameter  is  estimated,  from  the  figures  under  case  3  in  Tables  II  and  III,  to 
have  expanded  by  approximately  an  order  of  magnitude  by  the  time  AXu- 
Therefore,  in  this  case,  the  impulse  distribution  may  be  expected  to  be  a  fairly 
smooth  function  of  radial  distance  from  the  blast  center.  For  functions  of  this 
type  with  continuous  first  derivatives,  Kranzer  and  Keller  showed  that  only  one 
large  ripple  will  be  produced.  This  appears  to  be  in  agreement  with  the 
photograph  shown  in  Fig.  23. 

The  experimentally  measured  shock  peak  pressures  could  not  be  directly 
compared  with  the  theoretical  blast  model  predictions.  This  was  because  the 
rise  time  of  the  shock  probe,  at  1.5  ps  according  to  the  manufacturer's 
specifications,  was  not  quick  enough  for  the  predicted  time  constants  of  the 
shock  pulses,  as  shown  in  Table  IV.  For  this  reason,  the  shock  pulses  shown  in 
Figs.  16-18  are  probably  not  a  faithful  reproduction  of  the  actual  shock  pulse; 
Their  time  constants,  found  to  be  between  1.5  and  2  ps,  are  determined  mostly 
by  the  limitations  of  the  probe.  To  compare  these  results  with  the  model 
predictions,  it  is  necessary  to  allow  for  the  sluggishness  of  the  probe.  Assuming 
that  the  probe  behavior  may  be  modeled  as  a  linear  low  pass  filter,  the 
expected  peak  pressure  would  be  scaled  down  by  a  factor  approximately  equal 
to  the  ratio  of  the  time  constants.  The  experimentally  measured  peak  pressures 
are  compared  with  the  scaled  theoretical  predictions  in  Fig.  24.  It  is  seen  that  at 
0.03  MJ/m2,  the  intensity  of  the  laser  spot  is  less  than  the  minimum  required  for 
a  stable  blast;  this  is  supported  by  the  photographic  observations  which  show 
no  visible  surface  activity.  Nevertheless,  a  significant  level  of  acoustic  pressure 
was  produced.  The  mechanism  responsible  must  be  a  weak  blast,  i.e.,  one  that 
does  not  satisfy  the  maximum  entropy  condition.  The  results  at  1  MJ/m^  show 
good  agreement  with  the  equilibrium  blast  model;  photographic  observations 
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clearly  show  that  an  explosive  reaction  has  taken  place.  The  acoustic 
pressures  at  14  MJ/m2,  however,  are  considerably  greater  than  that  predicted 
by  the  equilibrium  blast  model.  Speculations  as  to  the  reason  for  this  excess 
pressure  will  be  given  in  a  later  section. 

Comparable  experimental  data  from  other  sources  have  not  been  found. 
There  have  been  a  large  number  of  published  reports  of  optoacoustic 
generation  experiments  in  the  interior  of  a  liquid  medium,  but  very  few  at  the 
surface.  The  most  extensive  data  set  is  that  of  Maccabee  and  Bell.''°-^^  They 
used  a  CO2  laser,  which  is  considerably  different  from  the  Nd:glass  laser  in 
terms  of  laser  wavelength.  Two  subsets  of  their  data  are  compared  with  the 
blast  model  predictions,  one  at  a  laser  spot  area  of  1  cm^  and  the  other  at  a  spot 
area  of  6  cm^,  as  shown  in  Figs.  25  and  26.  In  both  cases,  it  is  seen  that  the 
laser  intensities  used  in  the  experiment  were  too  low  for  an  equilibrium  blast. 
The  predicted  time  constants  in  the  equilibrium  blast  region  are  approximately 
2.5  |is  and  4  ps,  respectively.  The  measured  time  constants  were  in  the  region 
of  1  ps  for  both  cases.  It  is  deduced  that  the  mechanism  that  produced  the 
acoustic  pulses  must  have  been  a  weak  blast,  probably  due  to  a  number  of 
small  equilibrium  blasts  at  a  few  hot  spots  within  the  larger  laser  spot.  This 
would  account  for  the  reduced  time  constant.  This  is  also  consistent  with  the 
bright  flashes  observed  at  the  higher  intensity  levels. 

Comparison  of  the  predicted  and  measured  acoustic  energy  density 
spectra  were  made.  For  this  comparison,  there  was  no  need  for  any 
compensatory  adjustments  within  the  passband  of  the  acoustic  receiving 
system.  In  Fig.  27,  the  predicted  and  measured  spectra  at  1  MJ/m^  are  shown. 
It  is  seen  that  there  is  good  agreement.  Above  100  kHz,  the  experimental 
energy  density  spectrum  was  mostly  contained  within  the  leading  positive 
compression  pulse,  as  assumed  by  the  blast  model.  Below  100  kHz,  the 
experimental  spectrum  from  the  rarefaction  phase  of  the  shock  pulse  was  found 
to  be  higher  by  a  few  decibels.  In  Fig.  28,  a  similar  comparison  is  made  for  the 
14  MJ/m2  case.  Again,  it  is  seen  that  above  100  kHz,  the  experimental  energy 
density  spectrum  was  mostly  contained  within  the  leading  positive  compression 
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EXPERIMENTAL  RESULTS  OF  MACCABEE  AND  BELL  WITH 
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FIGURE  27 

COMPARISON  OF  THE  EXPERIMENTALLY  MEASURED  ACOUSTIC 
SPECTRUM  AND  THE  BLAST  MODEL  PREDICTION  AT  1  MJ/tn 
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FIGURE  28 

COMPARISON  OF  THE  EXPERIMENTALLY  MEASURED  ACOUSTIC 
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pulse.  It  is  also  seen  that  the  experimental  spectrum  level  is  far  higher  than  the 
equilibrium  blast  model  prediction. 


The  experimental  acoustic  spectra  from  various  observation  angles  q, 
where  the  angle  0  is  as  shown  in  Fig.  13,  are  compared  in  Figs.  29-31.  The 
results  for  the  0.03  MJ/m2  case  are  shown  in  Fig.  29.  The  spectra  appear  to 
decline  at  a  rate  of  about  20  dB  per  decade.  This  suggests  a  distributed  source 
function  rather  like  the  linear  thermoacoustic  source.  The  spectrum  in  the  0° 
direction  appears  to  be  the  strongest  except  below  150  kHz  where  the 
broadside  thermoacoustic  array  generates  a  significant  peak  at  near-horizontal 
angles  such  as  at  80°. 

In  the  1  MJ/m2  case,  as  shown  in  Fig.  30,  the  spectra  are  much  flatter, 
falling  off  at  a  rate  of  less  than  10  dB  per  decade,  suggesting  a  more  compact 
source.  Within  30°  of  the  downward  vertical  the  spectra  appear  to  be 
independent  of  angle.  At  85°,  there  is  a  peak  in  the  spectrum  at  around 
200  kHz,  indicative  of  a  weak  broadside  thermoacoustic  array. 

The  results  for  the  14  MJ/m2  case,  as  shown  in  Fig.  31,  are  most  puzzling. 
The  spectra,  at  0°  and  30°,  have  returned  to  a  slope  of  20  dB  per  decade,  while 
the  spectra  at  the  two  highest  angles  are  almost  flat.  This  is  indicative  of  a 
distributed  array  once  again,  but  probably  an  end-fire  array  rather  than  a 
broadside  array.  Further  speculation  will  be  given  in  a  later  section. 

It  was  also  found  that  at  1  and  14  MJ/m^,  there  was  often  not  just  one 
acoustic  pulse  but  a  train  of  several  pulses.  As  many  as  four  pulses  were 
detected,  as  shown  in  Fig.  32.  These  pulses  are  thought  to  be  bubble  pulses. 

D.  Optical  Results 

The  mechanism  for  initiating  the  blast  is,  however,  not  known.  For  CO2 
lasers,  with  which  many  studies  of  optoacoustic  generation  have  been  made, 
the  laser  wavelength  is  10.6|i.  At  this  wavelength,  the  absorption  coefficient  is 


0 


200  400  600  800 

FREQUENCY  -  kHz 


100C 


FIGURE  29 

MEASURED  ACOUSTIC  ENERGY  SPECTRA  AT  DIFFERENT 
OBSERVATION  ANGLES  AT  0.03  MJ/m^ 
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FIGURE  31 

MEASURED  ACOUSTIC  ENERGY  SPECTRA  AT  DIFFERENT 
OBSERVATION  ANGLES  AT  14  MJ/m^ 
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EXAMPLE  OF  OBSERVED  MULTIPLE  ACOUSTIC  PULSES 
GENERATED  BY  A  SINGLE  LASER  PULSE  AT  14  MjW 
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high  enough  that  simple  absorption  and  the  subsequent  increase  in  enthalpy 
will  adequately  account  for  the  blast  formation.  However,  for  a  Nd:glass  laser 
with  a  wavelength  of  1 .06p,  at  the  known  absorption  coefficient  in  water  of 
15  m'^  the  laser  surface  energy  density  required  to  cause  the  water  to  boil  is  in 
excess  of  20  MJ/m^;  therefore  simple  absorption  cannot  account  for  the  initiation 
of  the  blast.  It  has  been  postulated®  that  the  high  intensity  of  the  laser  radiation 
causes  the  water  to  undergo  a  dielectric  breakdown  which  results  in  a  sharp 
increase  in  the  absorption  coefficient.  The  increase  in  absorption  coefficient 
then  causes  the  energy  density  threshold  for  boiling  to  be  sharply  reduced  and 
the  remaining  laser  energy  is  then  absorbed  into  a  thin  surface  layer  which 
evaporates  and  starts  the  blast  process. 

In  an  attempt  to  observe  the  increase  in  absorption  coefficient,  an  optical 
probe  was  placed  in  the  water  at  a  depth  of  0.34  m  directly  under  the  blast,  as 
shown  in  Fig.  13.  The  laser  pulse  as  monitored  by  the  reference  probe  through 
the  beam  splitter  was  compared  with  the  observed  optical  pulse  at  the 
underwater  probe.  In  Fig.  33,  the  observed  optical  pulses  are  compared  for  the 
case  of  a  peak  intensity  of  0.1  TW/m^.  The  underwater  probe  was  at  a  depth  of 
33  cm.  Its  output  was  adjusted  to  compensate  for  the  absorption  loss  at  that 
depth.  It  is  seen  that  the  two  curves  are  very  similar,  and  so  it  is  deduced  that 
there  was  little  extra  attenuation  over  and  above  that  of  the  small  amplitude 
absorption.  At  a  higher  peak  intensity  of  0.4  TW/m^,  as  shown  in  Fig.  34,  there 
appears  to  be  significant  extra  attenuation,  and  the  underwater  pulse  appears 
to  have  lengthened,  possibly  due  to  secondary  emissions  at  other  wavelengths: 
the  PIN  silicon  photodiodes  are  sensitive  to  a  range  of  wavelengths  from  about 
Ip  to  0.4p,  with  a  peak  at  0.8p.  At  these  levels  of  intensity,  no  disturbance  was 
photographically  observable  at  the  surface. 

Higher  intensities  were  obtained  by  decreasing  the  laser  spot  size.  At 
intensity  levels  that  are  associated  with  the  equilibrium  blast  model  and  that 
produced  surface  disturbances  of  the  type  shown  in  Fig.  22,  the  secondary 
emissions  were  found  to  be  very  strong,  as  shown  in  Fig.  35.  At  similar  intensity 
levels,  but  where  the  surface  disturbance  corresponds  to  that  shown  in 
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FIGURE  33 

COMPARISON  OF  THE  INCIDENT  LASER  INTENSITY  WITH 
THE  OBSERVED  INTENSITY  BELOW  THE  SURFACE  AT  LOW  POWER 
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FIGURE  34 

COMPARISON  OF  THE  INCIDENT  LASER  INTENSITY  WITH 
THE  OBSERVED  INTENSITY  BELOW  THE  SURFACE  SHOWING 
INCREASED  ABSORPTION  AND  SECONDARY  EMISSIONS 
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FIGURE  35 

COMPARISON  OF  THE  INCIDENT  LASER  INTENSITY  WITH 
THE  OBSERVED  INTENSITY  BELOW  THE  SURFACE  SHOWING 
STRONG  SECONDARY  EMISSIONS 
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Fig.  23,  there  appears  to  be  a  significant  reduction  in  the  secondary  emission, 
as  shown  in  Fig.  36. 


From  the  above  optical  results,  it  is  seen  that  there  is  a  noticeable  increase 
in  absorption  as  the  laser  intensity  is  increased.  It  is  also  accompanied  by 
secondary  emissions.  However,  it  is  unlikely  that  the  cause  is  dielectric 
breakdown  because  the  absorption  appears  to  increase  gradually  with 
increasing  laser  intensity,  starting  at  0.4x10^^  W/m^.  At  the  highest  intensity 
used,  which  is  2x10’“^  W/m^,  the  electric  field  is  no  greater  than  2.7x10®  V/m. 
From  the  length^^  of  the  0-H  bond  (10'''°m)  and  its  energy^^  (2.6  eV),  the 
strength  of  the  H-0  bond  is  estimated  to  be  2.6x10''®  V/m.  Therefore,  dielectric 
breakdown  seems  to  be  an  unlikely  explanation. 

E.  Speculations 

The  foregoing  theoretical  analysis  has  produced  a  blast  model  for  the  laser 
reaction  at  the  air-water  interface  which  gives  plausible  results  that  seem  to  fit 
some  of  the  data.  Experimental  observations  of  the  optical  signal  have  shown 
that  there  is  some  extra  absorption  of  the  laser  at  high  intensities,  but  the 
breakdown  mechanism  that  must  be  present  to  initiate  the  blast  process  is  still 
unknown. 

It  has  been  observed  here  and  previously  by  others  that  sound  generation 
is  often  accompanied  by  a  strong  optical  emission  from  the  water.  Therefore,  it 
is  suspected  that  a  plasma  must  have  formed.  A  plasma  is  a  collection  of 
charge  carriers,  and  therefore  the  laser  radiation  must  be  causing  the  water  to 
release  charge  carriers.  It  is  known  that  the  ionization  constant  and  the 
conductivity  of  water  increase  with  temperature,  but  simple  absorption  of  laser 
energy  and  the  subsequent  increase  in  temperature  cannot  account  for  the 
charge  densities  required  for  plasma  formation.  With  CO2  lasers  at  a 
wavelength  of  10.6)i,  the  plasma  appears  to  form  in  the  ejected  vapor.  With 
Nd:glass  lasers  at  1.06p,  the  photographic  evidence  suggests  plasma  formation 
in  the  liquid.  Occasionally,  a  visible  plasma  trail  descending  into  the  liquid  is 
observed  photographically,  as  shown  in  Fig.  37. 
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FIGURE  36 

COMPARISON  OF  THE  INCIDENT  LASER  INTENSITY 
WITH  THE  OBSERVED  INTENSITY  BELOW  THE 
SURFACE  SHOWING  HIGH  ABSORPTION 
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FIGURE  37 

A  PLASMA  TRAIL  AND  IMPACT  CRATER  OBSERVED  AT  ENERGY 
DENSITY  14  MJ/m^,  STROBED  AT  0.8  ms  AFTER  LASER  PULSE 


The  energy  density  of  14  MJ/m  was  achieved  by  focusing  the  laser  beam 
into  an  area  of  only  0.023  mm^.  The  acoustic  output  obtained  was  found  to  be 
significantly  higher  than  that  predicted  by  the  equilibrium  blast  model.  Although 
the  blast  will  grow  rapidly  due  to  the  high  intensity,  the  blast  model  calculations 
indicate  that  it  should  not  be  as  effective  as  a  blast  of  equal  energy  that  starts 
from  a  larger  initial  area,  as  seen  in  the  comparison  between  cases  2  and  3  in 
Table  IV.  This  suggests  that  there  is  a  mechanism  that  is  quite  different  from  the 
hydrodynamic  relationships  on  which  the  equilibrium  blast  model  is  based. 
Another  significant  observation  is  that  this  unknown  mechanism  does  not 
produce  as  much  secondary  emission  as  the  equilibrium  blast.  This  suggests 
that  the  plasma  energy  has  been  diverted  away  from  the  emission  process. 
One  possible  explanation  is  that  the  plasma  may  be  coupled  to  a  sound  wave  in 
a  way  that  transforms  electronic  energy  into  acoustic  energy.^^ 

The  above  observations  suggests  two  interconnected  speculations.  The 
first  is  that  the  high  intensity  laser  beam,  perhaps  through  multiple  photon 
reactions,  is  capable  of  generating  a  plasma  in  water.  The  second  is  that,  at 
sufficiently  high  plasma  densities,  the  plasma  itself  may  act  as  an  acoustic 
amplifier,  channeling  electronic  energy  into  acoustic  energy.  The  shock  pulse 
driving  mechanism  would  be  electronic  rather  than  thermodynamic. 


IV.  LASER-WATER  INTERACTION 


A.  Electron  Charge  Carrier  Generation 

To  study  the  mechanism  for  plasma  formation,  a  simple  experiment  was 
conducted  to  measure  the  quantity  and  mobility  of  the  charge  carriers  that  the 
laser  radiation  produced,  if  any.  It  is  speculated  that  the  ionization  process  is  a 
gradual  buildup  process,  with  a  growth  rate  that  is  dependent  on,  among  other 
factors,  the  laser  intensity.  In  order  to  capture  the  process  in  its  earliest  stages 
of  development,  relatively  low  laser  intensities  were  used.  With  low  intensities, 
the  complication  of  thermodynamic  effects  may  also  be  avoided. 

The  apparatus  is  shown  in  Fig.  38.  Distilled  water  was  held  in  a 
cylindrical  test  cell.  The  cell  has  transparent  perspex  windows  to  allow  the  laser 
beam  to  enter  and  exit.  A  pair  of  cylindrical  electrodes  were  inserted  into  the 
water  to  measure  its  conductivity.  A  potential  of  approximately  6  V  was  applied 
between  the  electrodes  and  the  electric  current  time  history  was  recorded  on  a 
Nicolet  digital  oscilloscope.  The  laser  pulse  was  also  simultaneously  recorded. 
Due  to  the  large  electrical  capacitance  of  the  electrodes,  it  was  necessary  to 
use  a  fairly  long  laser  pulse.  The  laser  pulse  was  produced  by  spontaneous 
emission  from  the  Nd:glass  rod,  without  Q-switching.  This  produced  a  pulse 
lasting  around  500  ps.  The  laser  intensity,  at  40  MW/m^,  was  well  below  the 
known  breakdown  thresholds,  and  no  visible  reactions  were  observed.  The 
experiment  was  repeated  over  a  range  of  temperatures  from  20°C  to  80°C. 
From  the  time  histories  of  the  signals,  a  number  of  characteristics  may  be 
deduced,  including  the  yield,  the  mobility,  and  time  constant  of  the  carriers 
generated  by  the  laser  pulse.  The  results  are  summarized  as  follows. 

The  simplest  measurement  that  could  be  made  was  that  of  conductivity. 
The  conductivity  of  the  water  was  estimated  from  the  measured  dark  current,  the 
applied  voltage,  and  the  calculated  electrostatic  field  of  the  electrodes.  The 
conductivity  of  the  water  was  found  to  increase  with  temperature  within  the 
range  of  temperatures  considered.  The  increase  in  conductivity  per  unit 
increase  in  internal  energy  is  plotted  in  Fig.  39.  It  is  seen  that  the  thermal 
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FIGURE  39 

MEASURED  ENERGY  COEFFICIENTS  OF  CONDUCTIVITY 
AS  A  FUNCTION  OF  TEMPERATURE 
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coefficient  of  conductivity  is  approximately  constant  over  the  range  of 
temperatures  used.  The  values  are  comparable  to  estimates  based  on  the 
tabulated  values  of  ion  densities  and  mobilities  for  distilled  water.^^  These 
results  were  then  compared  with  the  increase  in  conductivity  due  to  absorption 
of  laser  energy,  measured  a  short  time  (60  ps)  after  the  start  of  the  laser  pulse. 
The  timing  of  the  measurement  was  important,  because  it  must  be  longer  than 
the  time  constant  (<1  ps)  of  the  electrical  circuit  as  determined  by  the 
capacitance  of  the  electrodes  and  the  parallel  resistance,  but  shorter  than  the 
decay  time  of  the  charge  carriers  produced.  The  laser  energy  absorbed  per  unit 
mass  was  estimated  from  the  measured  absorption  coefficient  of  the  laser  in 
water  at  the  intensities  and  temperatures  used.  As  shown  in  Fig.  39,  the  laser 
energy  produced  significantly  greater  increases  in  conductivity  than  simple 
heating  with  an  equal  amount  of  energy,  particularly  at  the  higher  temperatures. 

Having  found  that  the  laser  energy  produced  significantly  greater 
conductivity  increases  than  that  expected  from  simple  heating,  the  data  were 
analyzed  further  for  the  properties  of  the  active  charge  carriers.  It  is  known  that 
the  carrier  in  distilled  water  is  the  proton,^®  and  the  mobility  of  the  proton  is  well 
tabulated.  From  the  value  and  risetime  of  the  dark  current,  the  mobility  of  the 
proton  carriers  was  estimated  and  found  to  be  comparable  to  the  tabulated 
values,  as  shown  in  Fig.  40.  The  mobility  of  the  charge  carriers  generated  by 
the  laser  energy  was  also  estimated  and  compared  with  that  of  the  proton  in 
Fig.  40.  It  is  seen  that  the  mobilities  of  the  laser  generated  carriers  are  higher, 
by  a  factor  of  about  2000.  Since  electrons  are  about  2000  times  lighter  than  a 
proton,  it  is  guessed  that  the  laser  must  be  exciting  valence  electrons  to  a  state 
that  allows  them  to  act  as  charge  carriers.  They  are  probably  not  free  electrons 
as  in  the  case  of  a  metallic  conductor  but  more  akin  to  carriers  in  a 
semiconductor  that,  while  loosely  bound,  are  able  to  hop  from  one  molecule  to 
the  next.  Their  estimated  1/e  lifetimes  are  shown  in  Fig.  41 .  It  is  seen  that  they 
only  exist  for  a  short  time  about  100  ps.  Finally,  the  yield  in  terms  of  electrons 
per  photon  is  shown  in  Fig.  42.  It  is  seen  that  the  yield  is  fairly  constant  at  low 
temperatures,  but  it  increases  dramatically  with  temperature  above  55°C.  The 
mobility,  yield,  and  lifetimes  of  the  laser  generated  carriers  were  calculated  from 
the  measured  time  history  of  the  laser  generated  conduction  current,  particularly 


25  30  35  40  45  50  55  60  65  70  75 
TEMPERATURE  -°C 


—  TABULATED  PROTON  MOBILITY 
©  MEASURED  PROTON  MOBILITY 
♦  LASER  GENERATED  CARRIER 


FIGURE  40 

COMPARISON  OF  CARRIER  MOBILITIES 


the  initial  rate  of  current  increase,  and  the  conduction  currents  at  60  and  300  |is. 
They  are  related  through  a  set  of  nonlinear  simultaneous  equations.  The 
solutions  were  obtained  iteratively. 


B.  Plasma  Generation 

From  the  experimental  results,  it  appears  that  the  laser  photon,  at  a 
wavelength  of  1.06p,  is  able  to  liberate  electrons.  At  this  wavelength,  it  cannot 
directly  dislodge  an  electron,  and  therefore  it  must  be  energizing  the  electrons 
through  a  vibration  mode  of  the  water  molecule.  The  excitation  must  be  such 
that  electrons  are  energized  to  a  degree  that  they  can  significantly  behave  as 
charge  carriers.  Therefore,  it  is  speculated  that  a  plasma  can  be  formed.  The 
yields  within  the  measured  range  of  values,  however,  are  not  sufficient  to 
produce  a  plasma  of  a  sufficiently  high  density.  The  level  of  plasma  density  that 
is  considered  to  be  sufficient  will  be  defined  in  the  next  section. 

Given  that  the  photon  energy  is  efficiently  channeled  into  the  production  of 
electron  charge  carriers,  the  law  of  conservation  of  energy  dictates  that  the 
upperbound  of  the  electron  yield  must  be  governed  by  the  electron  energy 
increase  divided  by  the  available  photon  energy.  If  it  is  assumed  that  the 
electrons  are  not  entirely  free,  but  sufficiently  energized  that  they  can  hop  from 
one  molecule  to  the  next,  then  it  is  speculated  that  the  band  gap  energy  Eb  is 
the  difference  between  the  electron  affinities  of  the  H  and  OH  ions.  From  the 
electron  affinity  tables,^^  the  values  are  0.80  and  1 .83  eV.  Therefore,  the  value 
of  Eb  is  estimated  to  be  1.03  eV,  as  shown  in  Table  V.  At  1.06p,  the  photon 
energy  is  1.17  eV.  Therefore,  the  upperbound  of  the  electron  yield  is  estimated 
by  dividing  Eb  into  the  photon  energy,  giving  1.14  electrons  per  photon.  If  it  is 
assumed  that  the  observed  yield  trend  shown  in  Fig.  42  may  be  extrapolated  to 
higher  temperatures  until  the  upperbound  is  reached,  as  shown  in  Fig.  43,  then 
a  sufficiently  high  yield  would  be  achievable. 

Even  at  a  laser  energy  density  of  14  MJ/m^,  calculations  based  on  the 
small  amplitude  absorption  coefficient  indicate  that  the  heat  energy  absorbed 
would  only  raise  the  temperature  of  the  water  by  approximately  10°C  from  room 


TABLE  V 


EQUILIBRIUM  BLAST  MODEL  CASE  STUDIES:  DETONATION 


PLASMA  SOUND  SOURCE  PARAMETERS 

Laser  frequency 

Hz 

2.83E14 

Electron  affinity  of  the  OH  -ion 

eV 

1.83 

Electron  affinity  of  the  H  -ion 

eV 

0.80 

Difference 

eV 

1.03 

Photon  energy  at  l  .06|i 

eV 

1.170 

Uppertaound  electron  yield  per  photon 

1.136 

Critical  density 

m-3 

9.9E26 

Energy  density  threshold 

MJ/m2 

11 

CASE  STUDIES 

1 

2 

3 

Spot  diameter 

m 

0.0032 

0.00066 

0.00017 

Spot  area 

m2 

8.0E-6 

3.4E-7 

2.3E-8 

Incident  intensity 

W/m2 

5.6E-I-11 

1.3E+13 

2.0E-h14 

Energy  density 

MJ/m2 

0.03 

1.0 

14 

PLASMA  SOUND  SOURCE  CALCULATIONS 

UNIFORM  INTENSITY  MODEL 

Time  to  produce  plasma  cutoff  (tga 

)  PS 

19.1 

0.812 

0.054 

Time  for  complete  ionization  (tob ) 

ps 

197 

8.4 

0.55 

RANDOM  INTENSITY  MODEL 

Probability  of  avalanche  ionization 

% 

0 

100 

100 
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I 

I 
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temperature,  which  is  negligible  compared  to  the  temperatures  necessary  to 
produce  a  sufficiently  high  yield.  From  the  above  experimental  results,  it  may 
be  deduced  that  the  laser  photons,  at  a  wavelength  of  1.06p,  can  selectively 
excite  certain  modes  of  molecular  vibration  that  produce  conduction  electrons. 
Let  us  refer  to  these  modes  as  "conduction  modes".  If  a  temperature  can  be 
ascribed  to  the  degree  of  excitation  of  these  conduction  modes,  then  it  follows 
that  the  laser  photons  can  selectively  raise  their  temperature.  The  ratio 
between  the  temperature  rise  in  the  conduction  modes  and  the  overall 
temperature  rise  may  be  considered  as  a  "temperature  gain  factor".  Thus,  it  is 
postulated  that,  while  the  laser  energy  density  can  only  raise  the  overall 
temperature  of  the  water  by  a  few  degrees,  it  may  raise  the  temperature  of  the 
conduction  modes  by  the  equivalent  of  several  tens  or  hundreds  of  degrees. 
This  temperature  gain  is  expected  to  be  only  momentary;  its  duration  is 
expected  to  be  roughly  equal  to  the  lifetime  of  the  carriers  generated. 


In  order  to  obtain  an  explosive  reaction,  we  also  need  a  process  by  which 
a  plasma  can  concentrate  the  laser  energy  into  a  small  volume  to  generate  an 
intense  blast.  From  plasma  theory,  it  is  known  that  plasmas  have  an  e-m  wave 
cutoff  frequency  Wp,  called  the  plasma  frequency,29-36  given  by 

tOp  =  (Ne^/EomJ^^  ,  (80) 

where 

N  is  the  volumetric  carrier  density, 
e  is  the  carrier  charge, 

Eq  is  the  permittivity  of  the  medium,  and 

rrig  is  the  carrier  mass,  assumed  to  be  an  electron. 

All  e-m  waves  of  frequencies  lower  than  the  plasma  frequency  cannot 
propagate  through  the  plasma,  and  are  rapidly  absorbed.  The  critical  density 
No  of  an  electron  plasma  at  a  wavelength  of  1.06p  is  estimated  to  be 


9.9x1026  electrons/m3.  When  this  critical  plasma  density  is  achieved,  it  will 
cause  the  laser  absorption  to  abruptly  increase.  Subsequent  laser  energy 
arrivals  will  be  concentrated  in  a  thin  surface  layer  causing  local  superheating, 
and  hence  producing  the  gas  to  drive  a  blast.  The  plasma  will  also  radiate  over 
a  broad  spectrum  which  is  consistent  with  the  experimental  observations.  The 
time  toa  required  to  achieve  the  plasma  cutoff  is  given  by 

where 

cox  is  the  laser  frequency, 

Y  is  the  electron  yield  per  photon, 

Iq  is  the  laser  intensity, 

Ep  is  the  laser  photon  energy,  and 

ttoa  is  the  laser  absorption  coefficient  before  plasma  cutoff,  assumed  to  be 
equal  to  the  low  intensity  absorption  coefficient.  At  a  wavelength  of  1 .06|i,  aoa  is 
equal  to  ISm'''. 

Using  the  upperbound  electron  yield  of  1.14  electrons  per  photon,  the  time 
required  to  achieve  the  plasma  cutoff  for  the  three  numerical  case  studies  is 
given  in  Table  V.  After  plasma  cutoff  is  achieved,  further  laser  energy  arrivals 
will  be  concentrated  in  an  increasingly  thin  layer  of  the  medium.  The  absorption 
coefficient  ttob  is  now  given  in  terms  of  the  laser  frequency  co^  by^^ 

«ob  =  “X  [  Mo^o  •  (82) 

The  rate  of  increase  of  carrier  density  dN/dt  is  related  to  the  absorption 
coefficient  by 

dN/dt  =  Ya,blo/Ep  .  (83) 

Substituting  for  cOp  from  Eq.  (80)  into  (82)  and  for  ttob  from  (82)  into  (83), 
and  then  integrating,  an  expression  for  the  plasma  density  N  and  elapsed  time  t 
is  obtained. 


Setting  N  equal  to  the  number  of  available  conduction  electrons  Nee,  the 
additional  time  tob  required  to  achieve  total  ionization  may  be  estimated. 
Assuming  two  electrons  per  H2O  molecule,  the  available  density  is  calculated 
to  be  3.35x1 0^®.  As  shown  in  Table  V,  once  plasma  cutoff  has  been  achieved, 
it  is  only  a  matter  of  picoseconds  before  complete  ionization  is  achieved.  This 
process  is  called  avalanche  ionization. 

In  the  ionization  avalanche,  a  point  is  reached  where  the  energized 
electrons  outnumber  the  ground  state  electrons.  The  population  inversion  will 
cause  the  plasma  to  emit  e-m  energy.  Assuming  a  simple  two-level  energy 
state  system,  the  volumetric  source  density  Q  of  the  e-m  radiation  is  given  by®^ 

Q  =  (2N  -  Nce)EbAe  .  (85) 

where  tg  is  the  1/e  life  of  the  electrons  in  the  higher  energy  state. 

Multiplying  both  sides  of  Eq.  (79)  by  Iq,  we  obtain  Igtoa  on  the  left  which  is 
laser  energy  density  or  fluence,  and  the  terms  on  the  right  hand  side  are  all 
constants.  Thus,  the  above  theory  of  avalanche  ionization  predicts  a  definite 
avalanche  threshold  in  terms  of  the  laser  energy  density.  Assuming  the 
upperbound  yield,  the  threshold  is  computed  to  be  1 1  MJ/m^.  At  this  point,  it  is 
useful  to  make  comparisons  with  published  theories  and  experimental  results. 

It  has  been  suspected  by  a  number  of  researchers  that  free  electrons 
play  a  crucial  role  in  the  initiation  of  the  avalanche  process.  Bass  and 
coworkers, 3®-'^°  in  studying  the  breakdown  in  crystals,  concluded  that  the 
process  is  a  probabilistic  one  which  relies  on  the  existence  of  one  or  more 
"lucky"  electrons  that  are  able  to  resonate  at  the  frequency  of  the  laser  e-m 
waves,  and  thus  act  as  a  tuned  antenna  to  concentrate  laser  energy  into  a  small 
volume.  The  energy  concentration  then  spills  over  to  neighboring  electrons 
causing  them  to  also  resonate.  Thus  the  excitation  spreads  as  an  avalanche. 


There  are  several  published  measurements  of  the  optical  breakdown  of 
water, -^5  much  of  it  in  the  field  of  ophthalmic  laser  surgery.  In  most  studies,  it 
was  assumed  that  the  breakdown  threshold  is  one  of  intensity  or  electric  field 
strength.  In  each  case,  the  experiments  were  done  with  a  laser  pulse  of  a  fixed 
duration,  and  therefore  it  was  impossible  to  separate  intensity  (irradiance) 
dependence  from  energy  density  (fluence)  dependence.  The  most  recent  and 
thorough  study  is  probably  that  of  Docchio  et  al.^^  using  an  extremely  short 
laser  pulse  of  7  ns  and  a  range  of  spot  sizes;  the  spot  size  is  the  diameter  of  the 
narrowest  part  of  the  focused  laser  beam.  The  experiments  were  carried  out  at 
the  Center  for  Fast  Kinetics  Research,  The  University  of  Texas  at  Austin.  In  the 
absence  of  any  definite  theoretical  model,  they  took  the  precaution  of 
presenting  their  results  in  all  possible  forms  (intensity,  field  strength,  energy 
density,  etc.).  They  found  that  the  thresholds,  in  all  the  forms  used,  were  not 
constant  but  dependent  on  the  spot  diameter.  The  observed  intensity 
thresholds  from  several  sources  are  plotted  as  a  function  of  spot  diameter  as 
shown  in  Fig.  44.  The  estimated  dielectric  strength  of  the  0*H  bond  is  estimated 
to  be  2.6x10^°  V/m,  which  is  equivalent  to  an  intensity  of  1.3x10®TW/m^.  It  is 
also  noted  that  the  minimum  spot  diameter  is  approximately  equal  to  the 
wavelength  of  1.06p.  It  is  tempting  to  speculate  that  the  data  points  indicate  a 
trend  line  or  curve  that  terminates  at  the  point  of  intersection  of  the  smallest  spot 
diameter  and  the  molecular  dielectric  strength.  However,  the  trend  lines  in  the 
subgroups  are  not  consistent  with  this  view.  The  spread  of  data  points  is  also 
too  large  to  allow  a  regression  fit  of  any  significant  confidence  level. 

Plotting  the  breakdown  thresholds  in  terms  of  energy  density  appears  to  be 
a  small  improvement  in  that  the  spread  is  now  reduced  from  three  decades  to 
two,  as  shown  in  Fig.  45.  The  theoretical  energy  density  threshold  is  given  by 
the  product  of  the  laser  intensity  and  the  time  it  takes  to  achieve  plasma  cutoff, 
Igtoa.  which  was  computed  above  to  be  11  MJ/m^.  The  theoretical  threshold  is 
shown  in  Fig.  45  for  comparison.  Examining  the  experimental  data,  it  can  be 
seen  that  the  results  for  distilled  water  and  saline  solutions  are  almost 
indistinguishable.  The  data  for  tap  water  is  seen  to  be  very  variable,  probably 
due  to  the  wide  variation  in  the  quality  of  tap  water  and  particularly  the  presence 
of  dissolved  gases;  therefore  the  tap  water  data  are  best  ignored. 
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In  order  to  unravel  the  apparent  disorder  in  the  data  for  distilled  water  and 
saline  solutions,  let  us  speculate  on  the  possible  explanations  for  the  spot  size 
dependence.  The  theoretical  calculations  are  based  on  the  assumption  of  a 
uniform  laser  beam  profile.  This  is  clearly  not  the  case  in  practice.  Multimode 
laser  beams  have  a  more  or  less  random  beam  profile,  and  therefore  the 
energy  intensity  at  any  point  in  a  laser  beam  is  a  random  variable.  Focusing  the 
beam  increases  the  mean  energy  density  within  a  reduced  area  but  the  random 
nature  remains.  Therefore,  the  problem  may  be  treated  as  a  probabilistic  one. 

A  laser  beam  cross-section  may  be  divided  into  a  number  of  statistically 
independent  resolvable  areas.  Let  n  represent  the  number  of  independent 
areas.  The  random  variation  is  assumed  to  be  the  result  of  the  superposition  of 
several  modes.  Therefore,  for  large  n,  the  probability  density  function  of  the 
energy  density  at  any  point  may  be  approximated  by  an  exponential  density 
function.  For  a  mean  energy  density  Eq,  the  probability  p{Eo,n)  that  at  least  one 
resolvable  area  has  an  energy  density  greater  than  the  threshold  E,  is  given  by 

p(Eo,n)  =  1  -  (1  -  exp(-E,/Eo))n  .  (86) 

The  value  of  E,  is  11  MJ/m^  as  calculated  above.  The  size  of  the  minimum 
resolvable  region  is  determined  by  the  laser  wavelength.  Therefore,  at  any 
given  spot  diameter  do,  the  upper  bound  of  n  is  approximately  given  by 

n  {26JXf  .  (87) 

The  actual  value  of  n  is  not  easily  estimated  since  it  must  depend  on  a  number 
of  equipment  specific  factors,  such  as  aperture  size,  the  number  of  modes,  and 
the  lens  configuration.  The  data  collected  in  Fig.  45  suggest  that  there  is  a 
dependence  on  the  spot  diameter  do,  but  this  dependence  is  probably  different 
for  each  data  set.  Let  us  take  the  most  complete  data  set,  that  of  Docchio,  and 
see  if  a  plausible  relationship  between  n  and  do  can  be  found  that  is  consistent 
with  the  above  inequality.  Let  us  assume  that  n  may  be  approximately  related 
to  do  by  a  power  law. 
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(88) 


:  m>0  and  d^>do  , 

where  m  and  d,  are  constants  to  be  determined.  To  be  consistent  with  the 
above  inequality,  it  is  required  that  d^  be  greater  than  do  for  positive  values  of  m. 
Substituting  for  n  from  Eq.  (88)  into  (86).  an  expression  for  p(Eo,n)  is  obtained  in 
terms  of  Eq  and  do,  and  m  and  d^. 

An  attempt  was  made  to  fit  this  model  to  the  results  of  Docchio.  The 
results,  as  shown  in  Fig.  46,  show  the  energy  density  thresholds  at  0%,  50%, 
and  100%  as  measured  by  Docchio  with  distilled  water  and  saline  solution;  the 
0%  and  100%  thresholds  are  theoretically  impossible  to  measure,  and  therefore 
they  are  interpreted  as  <10%  and  >90%,  respectively.  A  reasonably  good  fit 
was  obtained  with  m=2,  di=0.6  mm  if,  in  addition,  the  energy  density  Eq  is 
assumed  to  have  a  small  offset  of  -0.6  MJ/m^.  The  offset  may  be  due  to  low 
intensity  leakage  of  laser  energy  before  and  after  the  Q-switched  pulse,  a  well 
known  problem  with  the  Q-switched  laser.  The  10%,  50%,  and  90%  thresholds 
obtained  from  Eq.  (86)  at  the  stated  parameter  values  are  compared  with  the 
measured  values.  It  is  seen  that  the  50%  theoretical  curve  could  be  fit  quite  well 
to  the  corresponding  measured  data  points,  while  the  data  points  at  >90%  and 
<10%  fall  above  and  below  the  90%  and  10%  curves,  respectively.  Therefore,  it 
is  seen  that  the  theoretical  model  can  be  fitted  to  Docchio's  data  quite  well. 

Let  us  now  apply  these  results  to  our  three  numerical  case  studies,  as 
shown  in  Table  V.  Using  Eqs.  (86)  and  (88)  and  the  above  values  for  the 
constants  m  and  d^,  the  probabilities  of  avalanche  ionization  for  the  numerical 
case  studies  were  computed.  The  results  were  found  to  be  0%,  100%,  and 
100%,  to  at  least  the  seventh  decimal  place,  for  cases  1,  2,  and  3,  respectively. 
These  results  are  consistent  with  the  experimental  observations,  since 
breakdown  was  always  obtained  at  energy  densities  similar  to  those  of  cases  2 
and  3,  but  never  in  case  1. 
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In  the  light  of  the  above  model,  it  appears  that  avalanche  ionization  may  be 
achieved  in  small  localized  areas  even  when  the  mean  laser  energy  density  is 
less  than  the  threshold  value  of  1 1  MJ/m^.  These  "hot"  areas  will  form  the  nuclei 
from  which  ionization  may  spread.  Once  the  cutoff  plasma  density  is  achieved 
in  one  small  area,  it  would  only  take  a  few  picoseconds  for  the  available 
electrons  in  that  region  to  achieve  a  population  inversion  and  radiate  energy. 
Thus,  the  plasma  will  redirect  the  absorbed  laser  energy  to  the  surrounding 
molecules.  The  increased  radiation  will  then  energize  the  electrons  in  the 
surrounding  molecules.  Thus,  the  cutoff  plasma  will  grow  outwards  from  the 
nucleation  sites  and  spread  over  the  laser  illuminated  water  surface.  This  will 
cause  a  big  increase  in  the  laser  absorption.  Subsequent  laser  energy  arrivals 
will  be  absorbed  and  concentrated  in  a  very  thin  surface  layer  causing  it  to 
expand  explosively.  This  forms  the  one-dimensional  stage  of  the  equilibrium 
blast  model.  Subsequent  reactions  can  be  expected  to  follow  the  blast  model 
developed  in  the  preceding  sections. 

This  model  provides  a  more  plausible  mechanism  for  triggering  the  blast 
process  than  dielectric  breakdown  of  the  water  molecule.  It  also  appears  to  fit 
the  observations.  As  shown  in  the  photograph  of  Fig.  22,  and  in  the 
corresponding  recorded  optical  emission  signal  in  Fig.  35,  there  is  evidence  of 
strong  optical  emissions  indicative  of  plasma  formation. 

E.  Plasma  Sound  Source:  A  New  Mechanism 

In  cases  where  the  mean  laser  energy  density  is  greater  than  the  threshold 
of  1 1  MJ/m^,  such  as  in  case  3,  avalanche  ionization  will  be  simultaneous.  The 
sequence  of  events  is  expected  to  be  as  follows.  In  the  initial  plasma  buildup 
period  (t  <  toa),  laser  energy  will  be  absorbed  by  the  water  column  within  the 
laser  beam,  down  to  a  depth  determined  by  the  absorption  coefficient.  At  the 
end  of  the  buildup  period  (t  =  toa).  a  thin  surface  layer  will  achieve  plasma  cutoff. 
A  few  picoseconds  later  ( toa<  t  <  toa  +  tob)  the  plasma  in  this  layer  will  achieve 
population  inversion  and  begin  to  radiate  energy.  Subsequent  laser  energy 


arrivals  will  be  absorbed  and  concentrated  in  this  layer  causing  an  explosive 
reaction.  The  shock  wave  produced  is  estimated  to  be  strong  enough  to 
significantly  compress  the  adjacent  water  layer  beneath  and  raise  its  plasma 
density  to  the  cutoff  level  by  virtue  of  the  reduction  in  volume.  With  the 
reradiated  energy  from  the  first  ionized  layer,  the  adjacent  layer  will  also  go 
through  a  similar  process  of  population  inversion,  energy  reradiation,  and 
explosive  expansion.  The  explosive  expansion  will  add  energy  to  the  shock 
front.  Therefore,  the  shock  wave  is  amplified  as  it  propagates  vertically 
downwards  through  the  energized  medium.  Similar  reactions,  in  which  a 
population  inversion  wave  is  preceded  by  a  shock  front,  are  known  to  form  the 
basic  mechanism  of  certain  types  of  lasers.'^®  The  effect  is  also  similar  to  that  of 
a  detonation  wave  traveling  vertically  downwards  through  a  shaped  charge; 
here,  the  energized  water  column  functions  as  the  shaped  charge.  The 
acoustic  effect  would  be  similar  to  that  of  an  endfire  array  directed  vertically 
downward.  It  is  speculated  that  this  mechanism  will  adequately  account  for  the 
extraordinary  acoustic  output  measured,  as  shown  in  Fig.  24,  at  the  laser 
energy  density  of  14  MJ/m^.  It  is  probable  that  the  process  may  be  modeled  as 
an  equilibrium  blast  in  which  the  energy  transfer  mechanism  is  not  thermal 
conduction,  but  secondary  emission.  The  secondary  emission  is  likely  to  be 
coherent  with  the  incident  laser  energy,  and  therefore  the  underwater  front  will 
be  guided  along  the  laser  beam.  This  appears  to  be  a  new  mechanism  of 
optoacoustic  sound  generation. 

There  are  a  number  of  observations  that  tend  to  support  the  new 
mechanism  hypothesized  here.  First  of  all,  the  photoconductivity 
measurements,  at  the  very  minimum,  indicate  that  laser  excitation  of  water 
produces  more  charge  carriers  than  that  expected  from  simple  heating,  and  that 
the  carriers  are  quite  different  from  the  proton  carriers  normally  found  in  water. 
The  results  suggest  the  possibility  of  electron  plasma  generation  through  a 
buildup  process.  Therefore,  electron  plasma  interactions,  rather  than  dielectric 
breakdown,  are  the  suggested  mechanism  for  initiating  the  explosive  reaction  at 
the  water  surface;  plasma  formation  here  is  taken  to  mean  the  liberation  of 
electrons,  while  dielectric  breakdown  is  the  breaking  of  the  0-H  bonds. 


Observations  of  the  surface  ripples  show  a  plume  at  14  as  in 

Fig.  23,  which  is  characteristic  of  a  strong  downward  impulse;  this  tends  to 
support  the  hypothesis  of  a  vertically  directed  detonation  wave.  This  contrasts 
with  observations  at  1  MJ/m^  which  do  not  show  any  plume,  as  in  Fig.  22. 

At  14  MJ/m^,  a  few  photographs  of  the  surface  show  a  luminescent  trail 
indicative  of  a  narrow  plasma  column  penetrating  into  the  water,  as  shown  in 
Fig.  37.  It  is  not  possible  to  determine  from  this  photograph  if  the  plasma 
column  ignited  spontaneously  at  all  points  or  propagated  into  the  water  behind 
a  shock  front  as  postulated  above.  The  acoustic  signal  directivity  may  be  used 
to  help  resolve  ihis  ambiguity.  Comparing  the  acoustic  signal  spectra  from 
various  angles,  as  shown  in  Fig.  31 ,  it  is  seen  that  the  spectrum  is  strongest  in 
the  vertical  direction,  particularly  around  a  frequency  of  100  kHz.  This  suggests 
a  source  that  behaves  as  a  vertical  endfire  array  whose  length  is  comparable  to 
the  acoustic  wavelength  at  100  kHz.  If,  in  contrast,  the  plasma  ignited 
spontaneously  everywhere  in  the  column,  it  would  have  the  characteristics  of  a 
vertical  broadside  line  array;  this  model  is  mathematically  similar  to  the 
thermoacoustic  model.  Calculations  for  the  vertical  broadside  line  array  show 
the  acoustic  energy  spectra  peaking  in  a  near-horizontal  direction,  as  shown  in 
Fig.  47,  which  is  quite  the  opposite  of  the  observed  directivity  as  shown  in 
Fig.  31.  Thus,  the  photoconduction,  photographic,  and  acoustic  data  all  appear 
to  be  consistent  with  the  new  mechanism  hypothesized. 

The  theory  of  Kranzer  and  Keller^^  was  used  to  obtain  a  relationship 
between  the  ripple  structure  and  the  vertical  impact  of  the  laser  blast.  An  optical 
model  of  the  photographic  imaging  system  was  constructed  to  relate  the  ripple 
structure  to  the  observed  grid  deformations.  Thus,  the  photographic  image  may 
be  inverted  to  give  an  estimate  of  the  magnitude  of  the  vertical  impulse  and  the 
depth  of  the  central  impact  crater.  Using  the  wave  theory  of  Kranzer  and  Keller 
to  fit  the  grid  deformations  observed  in  Fig.  37,  the  impact  crater  was  found  to 
be  long  and  narrow.  The  optical  model  reconstruction  of  the  grid  image  and 
impact  crater  is  shown  in  Fig.  48.  This  result  is  also  consistent  with  the 
hypothesis  of  a  vertical  plasma  column  penetrating  into  the  liquid  behind  a 
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shock  front.  From  the  photograph,  it  appears  that  the  penetration  of  the  plasma 
column  is  significantly  deeper  than  the  depth  of  the  theoretically  estimated 
impact  crater. 

The  new  mechanism  is  different  from  existing  theories  because  it  requires 
neither  a  direct  dielectric  breakdown  of  the  medium  nor  vigorous  evaporation 
through  the  direct  conversion  of  laser  energy  into  heat  to  account  for  the  sudden 
generation  and  release  of  mechanical  energy  from  within  a  small  volume. 
Although  plasmas  have  been  observed  in  connection  with  the  laser  generation 
of  sound  in  water,  they  have  always  been  considered  as  a  by-product  of  the 
superheating  of  water.  Here,  it  is  postulated  that  the  plasma  is  an  essential  part 
of  the  mechanism.  This  new  mechanism  does  not  replace  the  blast  model.  It  is 
in  fact  a  new  variation  of  the  equilibrium  blast  model  in  which  heat  transfer  is  not 
conduction  limited  but  secondary  emission  limited.  The  conduction  limited 
model  is  still  expected  to  be  valid  in  many  cases,  particularly  in  the  case  of 
optoacoustic  generation  using  CO2  lasers  where  the  absorption  coefficient  is 
high  enough  that  direct  boiling  of  the  water  will  be  the  dominant  process. 

A  new  mechanism  of  optoacoustic  generation  of  sound  is  put  forward.  It  is 
only  a  hypothesis  at  this  stage.  It  is  based  on  a  number  of  postulates  that 
remain  to  be  proven.  All  the  results  obtained  appear  to  support  the  hypothesis, 
but  more  direct  and  rigorous  experimental  proofs  are  needed. 
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V.  CONCLUSIONS 


The  objective  of  this  study  has  been  to  investigate  nonlinear  optoacoustic 
generation  processes.  The  initiative  for  this  project  arose  out  of  the  findings 
from  a  previous  ONR  funded  project  on  the  thermoacoustic  sound  source.  While 
it  was  found  that  the  thermoacoustic  sound  source  was  excessively  inefficient, 
the  nonlinear  optoacoustic  sound  source  appeared  to  have  some  potential  as  a 
noncontact  underwater  acoustic  projector.  The  nonlinear  processes  are  a 
collection  of  processes  by  which  a  high  intensity  laser  pulse  is  explosively 
converted  into  sound  waves.  Due  to  their  explosive  nature,  many  aspects  of 
blast  theory  are  applicable. 


Work  had  been  initiated  under  the  previous  contract  to  model  the  blast 
process  at  the  air-water  interface.  The  validity  of  this  initial  model  was  limited  to 
relatively  low  laser  energy  densities  where  the  acoustic  propagation  could  be 
assumed  to  be  linear.  In  the  present  project,  a  new  model  was  developed 
which  could  be  coupled  to  nonlinear  acoustic  propagation  theories.  The  new 
model  is  called  the  equilibrium  blast  model.  As  the  name  implies,  it  is  based  on 
the  hydrodynamic  equations  of  equilibrium.  The  conditions  at  the  shock  fronts 
are  determined  by  the  maximum  entropy  condition.  The  distribution  of  energy 
within  the  blast  is  assumed  to  be  through  thermal  conduction.  Experiments 
were  carried  out  to  test  the  model,  and  the  model  was  also  compared  with 
published  experimental  data  by  Maccabee  and  Bell.  Comparison  of  the 
theoretical  and  experimental  data  showed  that  more  than  one  process  was 
involved.  So  far,  three  distinct  blast  processes  have  been  identified,  the  weak 
blast,  the  thermal  conduction  limited  equilibrium  blast,  and  the  emission  limited 
equilibrium  blast.  The  model  developed  in  this  project  pertains  to  the  second 
process.  The  first  process  is  thought  to  be  a  corruption  of  the  equilibrium  blast 
process  in  cases  where  there  is  not  enough  intensity  to  fully  develop  the 
equilibrium  blast.  The  third  is  a  new  process  that  had  not  been  anticipated;  it 
was  discovered  by  experiment.  The  processes  occupy  mutually  exclusive 
regimes  in  intensity-energy  density  space,  as  mapped  in  Fig.  49.  The 
boundaries  shown  are  estimates  based  on  theoretical  analyses  and 
experimental  results  obtained  at  the  locations  indicated  on  the  map.  The 
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boundary  between  the  weak  blast  and  the  equilibrium  blast  is  a  rather  fuzzy  one 
that  is  determined  by  the  laser  spot  size,  and  possibly  certain  other  equipment 
specific  factors.  The  diagram  also  shows  that  the  results  obtained  so  far  only 
represent  a  single  cut  through  the  various  regimes. 

It  was  deduced  from  the  experimental  results  that  an  essential  ingredient  of 
the  equilibrium  blast  processes  is  plasma  formation  and  subsequent  avalanche 
ionization.  To  investigate  the  process  of  plasma  formation,  a  simple  experiment 
was  carried  out  to  study  the  ionization  of  water  in  the  presence  of  laser 
illumination.  Extrapolation  of  the  results  indicates  that  plasma  generation  and 
avalanche  ionization  through  a  build-up  process  are  plausible.  Through  the 
examination  of  laser  induced  breakdown  data  from  several  sources,  a  semi- 
empirical  model  of  the  probability  of  avalanche  ionization  was  formulated.  It 
was  concluded  that  electron  plasma  interactions,  rather  than  dielectric 
breakdown,  is  the  more  likely  mechanism  for  initiating  the  explosive  reaction  at 
the  water  surface;  plasma  formation  here  refers  to  the  liberation  of  electrons, 
while  dielectric  breakdown  is  normally  taken  to  mean  the  breaking  of  the 
molecular  bonds,  in  this  case  0-H  bonds. 

The  emission  limited  equilibrium  blast  regime  was  not  anticipated  but  was 
discovered  experimentally.  It  is  a  regime  that  is  deserving  of  further  study  from 
a  basic  physics  point  of  view  and  from  an  applications  standpoint.  From  a  basic 
physics  point  of  view,  it  is  interesting  because  it  probably  involves  an  electronic 
process  of  shock  wave  generation  and  amplification,  rather  than  the  better 
known  hydro-  or  thermodynamic  processes.  The  acoustic  shock  wave  appears 
to  draw  its  energy  from  a  plasma  column.  From  the  applications  standpoint,  it 
appears  to  be  the  most  efficient  of  all  four  optoacoustic  generation  processes 
identified.  Furthermore,  it  has  the  potential  to  be  tuned.  The  length  of  the 
plasma  column  is  expected  to  control  the  time  constant  or  resonance  frequency 
of  the  shock  pulse.  Thus,  the  frequency  content  of  the  acoustic  output  may  be 
controlled  through  the  length  of  the  plasma  column.  In  practice,  this  may  be 
achieved  through  the  proper  selection  of  a  laser  wavelength  with  the 
appropriate  absorption  coefficient. 


Our  experimental  results  of  the  emission  limited  equilibrium  blast  show 
acoustic  spectrum  levels  over  50  dB  greater  than  that  predicted  by  the 
thermoacoustic  model  and,  in  the  region  of  the  spectral  peak  around  50  kHz, 
30  dB  higher  than  that  predicted  by  the  conduction  limited  blast  model.  With  a 
laser  energy  of  only  0.7  J,  the  experimental  results  showed  a  measured 
acoustic  source  level  of  200  dB  re  1  pPa  at  1  m,  in  the  band  from  10  to  200  kHz. 
While  the  source  level  is  similar  to  results  obtained  by  Maccabee  and  Bell  with 
a  CO2  laser  of  a  larger  beamwidth,  the  CO2  laser  data  were  found  to  have  a 
spectral  peak  at  a  higher  frequency  of  185  kHz.  A  spectral  peak  at  50  kHz 
makes  this  mechanism  attractive  as  a  potential  sonar  source.  The  prospect  of 
being  able  to  control  the  frequency  through  the  length  of  the  plasma  column 
and  perhaps  achieving  even  lower  frequencies  is  worth  further  investigation. 
For  these  reasons,  the  emission  limited  blast  is  potentially  able  to  deliver  a 
breakthrough  in  optoacoustic  sound  generation  for  sonar  applications. 
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The  steps  used  in  processing  the  data  from  the  experiment  are  outlined  in 
Fig.  50.  All  programs  were  written  in  FORTRAN  and  executed  on  a  CYBER  830 
computer.  Experimental  data  were  collected  by  a  Nicolet  4094A  oscilloscope 
and  stored  on  disk.  The  data  interface  program  allowed  for  the  transfer  of  the 
data  from  the  Nicolet  to  the  CYBER  in  an  integer  format.  This  file  was  then 
normalized  through  the  use  of  a  normalization  program  and  output  in  a  floating 
point  format.  The  floating  point  data  are  then  processed  using  a  compensation 
program.  The  program  compensates  the  data  by  removing  the  transfer  function 
of  the  recording  equipment.  The  gain  and  phase  response  of  the  equipment 
was  measured  as  a  function  of  frequency  using  the  Nicolet  oscilloscope  and  a 
variable  oscillator.  Amplifier  and  filter  effects  as  well  as  hydrophone 
sensitivities  were  included  in  the  response  function.  The  output  of  the 
compensation  program  is  the  estimated  time  history  of  the  acoustic  pressure. 
Finally,  the  acoustic  spectrum  estimation  program  produces  a  spectral  density 
of  the  acoustic  energy  density.  Both  the  pressure  signal  and  the  spectrum  may 
be  plotted. 
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FIGURE  50 
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Photographic  images  of  the  water  surface  deformation  were  recorded.  A 
diagram  of  the  photographic  apparatus  is  shown  in  Fig.  14.  The  tubing  was 
designed  to  position  the  upper  and  lower  lenses  for  the  laser  beam  as  well  as  to 
align  the  camera  and  the  grid  imaging  lenses.  A  stroboscope  illuminated  a  grid 
pattern  which  was  projected  onto  the  water  surface.  A  timing  circuit  was  used  to 
control  the  time  delay  between  the  laser  pulse  and  the  flash  pulse  of  the 
stroboscope.  The  photographs  obtained  show  a  deformation  of  the  grid  pattern, 
which  indirectly  reflects  the  deformation  of  the  water  surface.  A  modeling 
process  was  undertaken  to  gain  a  better  understanding  of  the  photographic 
images  observed.  The  purpose  of  the  model  was  to  numerically  simulate  the 
photographic  image  so  that  a  relationship  might  be  established  between  the 
photographic  image  and  the  actual  surface  deformation. 

Light  rays  originating  from  points  on  the  grid  were  traced  through  the  first 
imaging  lens  to  their  image  points  on  the  water  surface,  reflected  from  the  water 
surface,  and  finally  traced  through  the  camera  lens  onto  the  final  image  on  the 
film.  The  origin  of  the  coordinate  system  was  fixed  at  the  intersection  of  the  grid 
and  the  ray  that  passes  through  the  center  of  the  first  imaging  lens  to  the 
intersection  of  the  laser  beam  axis  and  the  undisturbed  water  surface.  Basic 
vector  calculus  methods  were  used  to  trace  the  incident  and  reflected  rays.  The 
apertures  of  the  lens  determined  the  visibility  of  each  ray. 

The  equations  of  Kranzer  and  Keller32  were  used  to  model  the  surface 
disturbance.  Cylindrical  symmetry  was  assumed  about  the  axis  of  the  laser 
beam.  In  the  case  of  the  results  shown  in  Fig.  48,  the  process  was  modeled  as 
a  Gaussian  shaded  impulse.  The  magnitude  of  the  impulse  was  adjusted  to 
obtain  an  image  which  matched  the  grid  distortions  observed  in  the 
photographic  image.  The  result  was  a  surface  deformation,  at  a  time  delay  of 
0.8  ms,  consisting  of  a  bore  hole  0.26  cm  deep  surrounded  by  a  single  ripple 
with  a  peak  height  of  0.009  cm  at  a  radius  of  0. 11 5  cm.  The  disturbance  decays 
to  zero  at  a  radial  distance  of  0.675  cm.  The  profile  is  shown  in  Fig.  51 . 
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